IL 


NUOVO CIMENTO 


ORGANO DELLA SOCIETÀ ITALIANA DI FISICA 
SOTTO GLI AUSPICI DEL CONSIGLIO NAZIONALE DELLE RICERCHE 


Vou. XII, N. 3 Serie decima 1° Maggio 1959 


Diagonalization of Hamiltonian. 


L. M. GARRIDO and P. PASCUAL 


Junta de Hnergia Nuclear - Madrid 


(ricevuto il 26 Luglio 1958) 


Summary. — We present a general method to diagonalize the hamil- 
tonian of particles of arbitrary spin. In particular we study the cases 
of spin 0,3, 1 and see that for spin 3 our transformation agrees with 
Foldy’s and obtain the expression for different observables for particles 
of spin 0 and 1 in the new representation. 


1. — Introduction. 


FoLpy (!) gave a transformation that diagonalizes the Dirac particle hamil- 
tonian. In general, this diagonalization clarifies the physical interpretation of 
the mathematical formalism (2), as it has been shown by H. FESHBACH and 
F. VILLARS. CASE (*) extended this transformation to particles of spin 0 and 1 
using a hamiltonian formulation, whose covariance is not well manifested. 

We intend to present a transformation that does just that from a general 
viewpoint, and so, we start from the equations for particles of arbitrary spin 
as presented by UMEZAWA°s book Quantum Field Theory. Our transformation 
is valid for all spins and includes the precedent formulations. 


(1) L. L. Forpy and S. A. WoOUTHUYSEN: Phys. Rev., 78, 29 (1950). 
(2) H. FesaBAcx and F. ViLcars: Rev. Mod. Phys., 30, 24 (1958). 
(3) K. M. Case: Phys. Rev., 95, 1323 (1954). 
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2. — Unitary transformation. 


In the cases when the relativistic wave equation for particles of arbitrary 
spin . = Ù 
(B DI + m)y = 0 


may be written in the hamiltonian form 
Pop = (Py + MP)p (k = 1, 2, 3) 
where 


: i 
Po = — da; Dr = — 10; , eee A B= Ba; 


and $,, are the generators of infinitesimal Lorentz transformations defined as 


g°(B,b, "xi B,B,) ’ 


where g is a peculiar constant for each spin (4), we can find a general unitary 
transformation that diagonalizes the hamiltonian, 


HP, TMP: 


Such an expression for H has been studied by CASE (°). 
HEPNER (4) shows that there exists a unitary transformation 


S = exp [3ixg}] 
such that 
S28, = Rf, 


Sp Se Ba 


under which H becomes 
H'=S-HS=f,p,+mf,. 
Our problem will be solved if we evaluate the transformation R that makes 


R-H'R = Es; 
where E the total energy 
E=Vpim. 


(4) W. A. Hepner: Phys. Rev., 84, 744 (1951). 
(5) K. M. Case: Phys. Rev., 100, 1513 (1955). 
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To start with we will consider the simpler case when p, = p,; — 0; we call 
R, the corresponding transformation. 
The transformation R, is a rotation in the plane (1, 4) whose infinitesimal 
expression is 
OR, = 1+ Sud 
and, therefore 
fet eXp)| Si.) 


where © is the angle through which we have rotated 


D tbe. 


Pi. 


The infinitesimal rotation in the plane determined by the axis f, and an 
arbitrary unit vector n, is 


ÔR, =1+ 8, bo 


and, therefore, the finite rotation that we look for is 


ues 1 p 
i = exp [su À ee tg- zi, 


where 
B= 4/ D2. 


Reversing the transformation S the hamiltonian will be diagonalized by 


A= SES = exp|—ig BrPs tg! P aE 
p m 


that in the case of a Dirac particle g=+ yields Foldy’s transformation (1) 


L g uP to! 4 


Ay = ex : 
3 P 2 p m 


3. — Spin 0 and 1. 


Spins zero and one correspond to g=1. For a Duffin-Hemmer particle 
represented by the equation of motion 


(2,6, +m)y = 0, 


B uByPa + B,B,B, = By 46 ae à Bò vu” 
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there exists (55) a hamiltonian formulation as follows 


ot 


. oy 6 
H = ap + mb: , Hy =i, 7 
and the initial condition is | 


(HB, — m)y = 0 ’ 
where 


y, as Wa ? S = Pies ES BP 7 


Pr = — 0, , Op = — Sym (k, l, m cyclic). 


The hamiltonian is diagonalized by 


A =exp| (EP ig Zl, 


so that 
(I) AHA = EB,, E=Vp?+m, 


and the initial condition becomes 


siate tti 


{Ps(a-p), — mPi+m}p=0, 


where y’=A-1y. It is quite easy to see that it eliminates the zero energy 
states. 
Expression (I) may be expanded 


1 i 
RTE (8-p)?—z (BP), 


ù de dt = 
Ca IR EE a e TE O 


according to the formulae (A.1). 


4. — Dynamical variables. 


We call in the new representation O’ the operators corresponding to O- 
in the old one 
O = AOA 


and so, with the help of the formulae of the Appendix, we find Table I. 

As it has been shown by Foupy (1) for the Dirac particle the mean spin 
operator and mean orbital angular momentum are also in these case separately 
constants of motion. 


(5) N. KEMMER: Proc. Roy. Soc., A 173, 91 (1939). 
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TABLE I. 


Operator in the old 


Operator in the new 


representation representation 
È : i(B-p)+(E+m) i 
Position | x x'—x + EE mp (Bp)p- 
TEE ME 
F E(E+m) pas ne 
Momentum | p=—iV p'= p 
Hamiltonian | H=a-p+ mp, He = eh) By 
elocity t= waa +s, Pa EE+m) (ap) 
Orbital angular| (x xp) (x'Xp) = (axp)+ 
momentum | 1 
| ae Em P xo) XP aa (B Xp) 
Soi | 5 e 1 (on 
| oe? Am PP 
| 1 
| #07 (Bxp) 
| 
Projection | __ ap+mp, : 
operator | aoe E Ao Ps 
| i(Bp)-- (E+m) 
eye | ; PP} 3 5 He, 
Mean position |X = x+ EXE4m)? (B-p)p+| X'= x 
1 B 
| di E(E+m) (Sa E 
| 
| 3 
Mean velocity | X = n X' a 
Mean orbital x n: 1 Xx p' 
angular (A Xp) = (x xp) + E(E+m) (A'xp’) = (x xp) 
momentum 1 
(p Xo) xp — —(B xp) 
M i 1 
ean spin er + 
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The eigenvalues of the velocity operator «=a are +c and 0. This last 
one is eliminated by the initial condition. This instantaneous velocity is due, 
as in the case of Dirac particles, to a rapidly oscillating motion (Zitterbe- 
wegung) since by means of 0 =i[H, 0] we deduce_ 


a = Hp — E?a 
and 
ü —-— Fa, 


differential equation whose solution is 


a = Aexp[—i£Et] + Bexp[+i£t] ee 


where À and B are integration constants, expression whose interpretation is 
well known. 

If we divide the spin o’ into parts parallel and perpendicular to the mo- 
mentum by means of 


dl 
p° 


; Or. 


we get 
m 


J I 1 5 
T= Ot FFL (PB XP): 


It is easy to see that the expected value of (B Xp) is zero and therefc-e 


m 


/ 


TONDE 
pio 


La 


(o) = (o) + 


APPENDIX 


The formulae used for all these calculations are: 


(B .p)*” =p (Bp) (BP = p°(p-p) 

[x,(B-p)] =p [x, (B-p)?] = (8 -p)8 + B(B-p)} 

[B.(B-p)] =uUpxc)  [B, (B-p)?] = i{(B-p)(p xe) + (pxo)(B-p)} 
(A.1) | [(ap),(B-p)]= i8,p° [(ap),(B-p)?|= 2ip°(B-p)f.+ p*(a-p) 

[Ps,(B-p)] =—i(a-p) [fs, (B-p)?] = Bap? — 2(B-p)fa 

[s, (B-p)] =—i(Bxp) [o,(B-p)?] =—i{(B-pP)(BXxp)}+(BXp)B-P} 

[o, (B-p)] = ip fa, (B-p)?] = ip{(8-p)fa+ Bi(B-p)} ; 


[(B-p), (8 Xp)] = i{(o-p):p — op} , 
(8B -p)B(B-p) = (B-p)p, 
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(B-p)(ap) = — i(B-p)°f., (B-p)°(ap) = —ip°(B-p)f., 
(B-p)(Bxp)(B-p)=0, 
(B-p)(a-p)(B-p) =0, 
(8B -p)P.(8 -p) =0. 
and a representation of the Duffin-Kemmer algebra is obtained from those 


matrices given by, UMEZAWA (7) that we will call B, by means of the following 
unitary matrices: 


Spin 0 
1 1 
ae 0. 0 0 
V2 2 | 
rie OT 07 0) 
D 0 0 0 0|. 
0 DENZA dx 
: Lei 
V2 V2 | 
Bu = U*B,U 
ae a 6 (ag 
DI 14/2 
0 Jalla 0 
A NT 030 
i V2 ve 
<< agg ee ee 
Pa v ae | © 0 
Va V2 i pena x 
0 0 V2 V2 
0 0 0 
fe | 
/2 ini 0 
| 0 
. fo | 0 —1 
| | va | 
: à ISS Pc ER PR Rd 
i= 2 >) aaa 3 STRA [2 Pa | ? 
0 0 : | | 
| 
Li E 0 | 0 0 
i —i. | 
V2 V2: | 


(7) H. Umezawa: Quantum Field Theory (1956). 
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So 


0 


0 0 


0 
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1 
1 0 0520 
1 | 
Re is 
BEN = BA UR 
—1: 
0 0 0 0 
0 0 0 0 


RIASSUNTO (*) 


Presentiamo un metodo generale per la diagonalizzazione dell’hamiltoniana delle 
particelle di spin arbitrario. In particolare studiamo i casi degli spin 0, 3, 1 e vediamo 
che per lo spin + la nostra trasformazione si accorda con quella di Foldy ed otteniamo 
l’espressione, per differenti osservabili, per le particelle di spin 0 ed 1 nella nuova 
rappresentazione. 


(*) Traduzione a cura della Redazione. 
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On the Emission of Low Energy Photons 
from a Quantized System. 
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Istituto di Fisica dell Universita - Torino 
stituto Nazionale di Fisica Nucleare - Sezione di Torino 


(ricevuto il 4 Agosto 1958) 


Summary. — We solve here the problem of the emission of photons from 
a general quantized system when the reaction of the photons on the 
system is neglected. Till now this problem had been solved only in the 
case of the emission from a classical system. We deduce a general ex- 
pression for the average number of photons which are emitted when a 
time dependent perturbation produces a transition between two arbitrary 
states of the system. The discussion of this expression leads to the con- 
clusion that the average number of low energy photons emitted in any ex- 
periment which can really be done is finite: this conclusion generalizes 
a result previously obtained by the author in the case of a classical source 
to the case of quantized source. Moreover, neglecting the reaction of 
the photons on the charges, we show here that the average number of 
emitted photons is equal to the first term of the expansion in powers of e 
of the emission probability of one photon. The application of this con- 
clusion to the case of a transition between two definite momentum states 
shows that the Bloch-Nordsieck divergence and the infrared divergence 
of the perturbative calculation have their origin in the divergence of the 
same mathematical expression, although they have an entirely different 
physical meaning. Finally we show that the problem studied here in 
general has no meaning in the case of a time independent perturbation. 


1. — Introduction. 


The problem of the emission of photons of sufficiently low energy may be 
studied independently from the perturbative expansion using an approximation 
method in which the reaction of the electromagnetic field on the charges is 
neglected. 
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Such an approximation was introduced firstly by BLOCH and NORDSIECK (1). 

Till now, however, this method has been applied rigorously and extensi- 
vely only in the case in which the photons are emitted from a classical system 
of charges (?). For this classical case there are-several general and elegant 
treatments (3). The most simple and general of them is the one of THIRRING 
and TOUSCHEK (*4), which uses Heisenberg picture. 

The most known result concerns the case of a charge undergoing a transition 
between two states of motion of constant momentum; in this case the average 
number of emitted photons results infinite (Bloch-Nordsieck divergence), 
whereas the probability that a finite number of photons is emitted is found 
to be zero (1). 

In a previous work (4) I have studied the expression of the average number 
of photons emitted by a general classical system of charges, and I have shown 
that this number is always finite, provided the charges are confined in a finite 
region of space. Therefore I have concluded that in every experiment which 
can actually be done the number of emitted photons is finite, so that the 
divergence found by BLOCH and NORDSIECK has its source in the fact that 
the particular problem in which it appears corresponds to an experiment which 
cannot be made, because an infinite region of space is needed. 

The main purposes of the present work are: 


1) To give an exact meaning to the problem of the emission of low energy 
photons from a « given » quantized system, when the reaction of these photons 
is neglected. 


2) To give a general treatment of this problem, using Heisenberg picture. 


3) To apply the result to the case of the transition of a charge between 
two definite momentum states, and to show rigorously that the result found 
in the classical case holds also in the quantum mechanical one, provided the 
duration of the transition is neglected. 


4) To show that the number of photons emitted in an experiment which 
can really be performed is always finite, extending the result previously ob- 
tained in the classical case to the quantum mechanical case. 


(1) F. BLocH and A. Norpsinck: Phys. Rev., 52, 54 (1937). 

(2) Only a very particular case of emission from a quantummechanical system was 
treated rigorously by Frerz and PauLI (W. PAULI and M. Frerz: Nuovo Cimento, 15, 107 
(1938)). See also, J. M. JAUCH and F. RonRLICH: Helv. Phys. Acta, 27, 613 (1954). 

(3) W. THIRRING and B. TouscHEK: Phil. Mag., 42, 244 (1951); Roy J. GLAUBER: 
$ hys. Rev., 84, 395 (1951); J. SCHWINGER: Phys. Rev., 91, 728 (1953). 

(4) See also: R. AscoLi: Nuovo Cimento, 2, 413 (1955). 

( ) See also: Res Jost: Phys. Rev., 72, 815 (1947); R. AscoLi: Nuovo Cimento, 
2, 1 (1955). 
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5) To prove that a simple connexion exists between the results of this 
treatment and the result of the perturbative calculations. 


6) To discuss the connexion between the Bloch-Nordsieck divergence 
and the infrared divergence which arises in the perturbative calculation (6). 


2. — Considerations holding in both the cases of a classical and of a quantized 
source. 


Let us consider the equations of motion of an arbitrary system in Heisen- 
berg picture. 

The approximation treated here consists in the following proceeding; in 
the equations of motion of the charges we neglect the term of interaction with 
the photons of energy lower than some fixed value e in a given reference 
system, whereas nothing is neglected in the equation of motion of the electro- 
magnetic field (7). 

Let us call 


— Degrees of freedom e: all the degrees of freedom except those of the 
photons of energy lower than e; 


— Degrees of freedom p: those of the photons of energy lower than e 
(« low energy photons »). 


Then the problem is solved in two steps: 


1) The motion of the degrees of freedom e is studied as if the degrees 
of freedom p would not exist at all; that is the degrees of freedom e are treated 
as belonging to an independent system. We call it system e or « external 
system », because it corresponds to an external source in the classical case. 


(6) A preliminary account of some results of the present work has been published 
on Nature with a synthetic deduction. See R. AscoLi: Nature, 179, 727 (1957). 

(7) A simple formulation of this approximative method from the beginning is only 
possible using Heisenberg picture. In fact we are using an approximation in which 
the influence of the low energy photons on the charges is neglected, whereas the in- 
fluence of the charges on the low energy pbotons is taken into account. 

Therefore both the Schrédinger and the interaction pictures are not appropriate 
to deal with such a problem, because in the interaction Hamiltonian the influence of 
the field on the charges and the influence of the charges on the field cannot be sepa- 
rated. So, using these pictures it is not possible to face the problem clearly from the 
beginning, and it is necessary to make the required approximations at the appropriate 
points during calculation. 

On the contrary, in the Heisenberg equations of motion the influence of the field 
on the charges and the influence of the charges on the field are separated from the 
beginning. 
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In such a way an expression is obtained for the density of electric current j(x) 
as a function of only the quantities referring to the system e. 


2) The density of electric current j(x) so obtained is introduced in the 
equations of motion of the low energy photons and these equations are solved. 


We suppose here that the first part of the problem has been solved in some 
exact or approximate way. 

We observe that the expression of any operator O° referring to the degrees 
of freedom e at any time will contain only the quantities referring to system e, 
so that 0° will commute with any operator 0? and with any state vector |g} 
referring only to the low energy photons 


(1) [O- 0" — 0, 
(2) 0-19) —1p)0. 


The second part of the problem consists in the treatment of the equation (8): 


(3) OA=—j, 


where j is supposed to be known and where we are interested only in the com- 
ponents of the field A(x) of energy smaller than e. 
Let us define A(k) and j(k) through (°): 


af fi : % 
(4) AG) = Taj | dk exp [é(ka)]6(k2) A (k) , 
d Pe PL Dae ria 
(5) ie) =a | dk exp [i(har)]j(0) 


Defining in the usual way incoming and outgoing field operators A") 
and A(z), it is easy to show that the simple relation holds: 


(6) A™ (k) = A°(k) + de(k)j(k), 


where A™(k) and A°*(k) are defined analogously to A(k) by (4) and e(k) —1 
for ky > 0, .e(/) = — 1 for, ED (EU) 


(8) A and j denote fourvectors; e= 1; A= 1 in the following. 

(9) (kx) denotes fourdimensional scalar product. 

(19) See, for example, the work of THIRRING and TOUSCHEK (*) or the work of 
the author (4). 
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Introducing the quantization of the electromagnetic field, the observable 
to which we shall be mainly interested is the number N of photons, whose 
momentum % lies in a volume V, of the momentum space, and whose pola- 
rization vector is e(k). We have 


1 * 
(7) N= igi, [how A, (hy AL), 


Mi 
where A,(k) = (A(k)e(k)). 
Correspondingly we define the number N™ of outgoing photons using the 
outgoing field A° instead of A. 


We shall call |0> the vacuum state of the low energy free photons, |1ke> 
the state of one free photon of momentum % and polarization e. Then it is 


(8) A® (k) 


Os 0. 


3. — Precisation of the problem for the case of a quantized source. 


As it was remembered in the Introduction, till now only the case in which 
the external system is classical had been treated in general. 

Let us now suppose that the external system e is quantized. In this case 
it is necessary to specify what a kind of quantum-theoretical problem for the 
external system is considered. We shall first consider the case of a transition 
caused in the external system by a time-dependent perturbation of finite dura- 
tion between two eigenstates of an unperturbed energy. In this way an un- 
perturbed and a perturbed external system are defined. 


— Unperturbed external system. — We call x a complete set of observables 
referring to the external system commuting with each other and with the 
unperturbed energy, x, (r=1, 2,...) a set of eigenvalues of these observables, 
|x,» the corresponding eigenstates in Heisenberg picture, H, the energies of 
such eigenstates. We shall also introduce the corresponding eigenstates in 
Schrédinger picture |«,(#)>. Then it is: 


(9) [ar (a0)> =|a,> exp[— iE], <a (ao) | = exp [1H 9] <ct, |. 


— Perturbed external system. = We suppose that the origin of the time x, 
has been fixed so that the time-dependent perturbation acts during the finite 
time interval 0<a,< T. 


We define in the usual way the incoming and outgoing eigenstates |i"), 
|x®*) of the total energy in Heisenberg picture. Then we introduce the cor- 
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responding eigenstates in Schrédinger picture |" (x), |«2"*(a))> and the uni- 
tary transformation operators in interaction picture U*(x%,) and U° (x) (#4) re- 
ferring to the external system, so that: 


(10) : lon (ay) > = U2 (a) |; (%)> = > Us, (20) |a; (Zo)> » 


(11) las" (00)) = UE (20) |oer(@o)> = X UL, (00) |a(20)) - 


The S-matrix for the external system is then defined by 


(12) Si, = as" |a") 
and it is 
(13) U* (to) = U* (a)(S*)* . 
Let us choose the different pictures so that they coincide for 7 = — oo. 


Then for a perturbation acting only from the time 0 to the time 7 we have 
(14) S°= U°(F) = U2(0); 

fora =< 0 

UD): UT) = 15° US (ey) S045 Vi) (0 UO ee 
one di 

COREA RE VII (OTT UG (EU 


Finally we define incoming and outgoing states for the total system. We 
call |x,0'"> the state of the total system in which for 2, > — co the external 
system is in the state |x,), while there is no low energy photon present. Ana- 
logously we define |a, 0%, |a,1 kc°*). 


4. — Calculation of the average number of emitted photons. First part. 


The most interesting quantity is the average number N,; of photons emitted 
with a given polarization in a given volume of the % space when the external 
system has undergone a transition between a state |x,) and a state |«,>, sup- 
posing that initially no low energy photon is present (72). 

Then the state of the experiment is |a, 0”). 

An expression for N,, may be found as follows. 

(11) See for instance, B. A. LippMANN and J. SCHWINGER: Phys. Rev., 79, 469 (1950). 

(2) In Appendix I is shown that the corresponding probability Wy,, for the emis- 
sion of N photons may be immediately calculated, when N,, is known. 
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Let us introduce the observable de, having the value 1, when a measu- 
rement of the observables « at the end of the experiment gives the value &, 
(that is when at the end of the experiment the external system is in the state 
|a,>), the value 0 when the external system is in an eigenstate of the unper- 
turbed system different from |x,). 

Let us consider the expression 


(17) Car OM | NG out, |a, 0"). 


This is the expectation value of the observable N°"3,cu, in the state 
|x,0">» of the experiment. Let us now measure N and « at the end of the 
experiment. Then the observable N°"ò cut, has the value of N when the final 
state of the external system is |«,>, 0 when it is another state of the set 
|x;,>. Thus its expectation value (17) gives the average number of photons 
emitted in the transition r +s, multiplied by the probability that such a 
transition occurs. 

So N,, is obtained dividing the expression (17) by |S%,|?: 


C7 (DE | Wo" Ont. |a, Oe 
EAE 


(18) N,, — 


This expression now may be transformed in such manner as to eliminate 
from it all the quantities referring to the field of the low energy photons. 
Using (7) we obtain from (18): 


i 4! . 
e I | dk 3 (k*) <a, Où | ASE) AK) Ba, 


‘Si, 2a) AE 


Using (1) with O° = d,t,, O0? = A°"(k) and (6) we may transform the ex- 
pectation value Æ in (19) as it follows: 


(20) | B= <a 0"|A5"*(k)A°"(k)d 0, |æ, O™> = 


= (ar 0%] (AY (H) — aj" (6) dont, (AE) + ij.(k) |a, 0") . 


At this point we suppose to use the special Heisenberg picture which co- 
incides with the Schrédinger picture as 7 > — co. 

Wbeniaz0 = |\c,>|O>. 

Using (8) and (2) we obtain 


(21) E = (O|O><c%, | 72 (le) dove Je (k) ar) = Car | Fo (H) ou, Fel) [oe > 5 


13 - Ji Nuovo Cimento. 
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where we may return to a general Heisenberg picture in the last expression. 
So we have obtained an expression containing only quantities referring 
to tke external system. 
We have further: 


a 


(22) E=) 2 |]? (h) dont, ot) Cat |, (k je = 


2 


Co |4*( k) last. COSTA k) ) |a = | {ae (k) |a) | 
and, finally, substituting into (19), 


= 1 1 
23 Nes dk 5(k2) | Co" |j im 
"oa sae aa EE) esata 


2 


TEis is a general expression which permits to calculate N,, whenever the 
quantities referring to the external system are known. 

An important feature of this expression is that it is proportional to €, 
e being the coupling constant of the charges with the low energy photons (!). 


5. — Comparison with the result of the first order perturbative calculation of 
the probability of emission of one photon. 


The expression (23) may be used to calculate or to discuss directly the 
average number W,, of emitted photons. 

However we prefer firstly to compare it with the expression of the proba- 
bility W,,, of emission of one photon, when the transition >» —s occurs. 

This quantity is given by 


Al 
LS. eo 


(24) Wa | dk 6(k?) | <a, 1ke™ |x, O |: 


The matrix element M in the expression (24) can be transformed as it 
fellows: 


(25) Me Cad he a, OS = <a, "| Al" (he) |e, ON = Ka O°" 9s) eee 


where j,(k) = (j(k)e(k)) and we have used (6) and (8). 

Now let us call e the coupling constant of the charges with the low energy 
photons. Then let us suppose to calculate W,,, with the perturbative method, 
using e as expansion parameter, and to retain only the first term W°” of the 


Irs 


(>) This property has been already pointed out and used in my paper (5). 
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expansion. ‘The result of such a calculation will be obtained expanding the 
matrix element M which appears in (24) in powers of e, and retaining only 
the first term M of the expansion. 

Let us first expand |x,0") in powers of e. The first term of the expan- 
sion is obtained putting e= 0, that is eliminating the interaction. Now with- 
out interaction the field of the low energy photons at any time is still in the 
vacuum state, while the external system is in the state |«™) in any case, so 
that we obtain: 


1x, O7» = |a) |O) + terms of higher order in e. 


Then we use the Heisenberg picture which coincides with the Schrédinger 
picture at the end of the experiment. Then <«,0™|—= <«°"|<O|. Now ob- 
serving that <x°"*|< O| does not depend on e, and that j,(k) is simply propor- 
tional to e, we have, using (2): 


(26) MOKONA). 
At this point we may return to a general Heisenberg picture, since the last 
expression is invariant for transformation of picture. 


Substituting into (24) we obtain: 


1 1 
(27) Wi rs fa 6 (k?) | os" 


[| Ss-|? 27)? EA E 


The comparison of this expression with (20) leads to the important result (14): 
(28) Ne Mas ’ 


that is, when the reaction of the photons is neglected, the average number 
of photons emitted in a transition is equal to the first term of the expansion 
in powers of e of the emission probability of one photon in the same transition. 

The consequences of this result will be discussed later. Here I observe only 
that this result reduces the calculation of N,, to a perturbation problem. 
However such a perturbation problem till now has been treated only in par- 
ticular cases, neglecting the duration of the transition. 

Here a general rigorous method of calculating expression (23) will be given, 
with the purpose of discussing the result in the most interesting cases. 


(4) This result has been reported in (*) where it was deduced without direct 
calculations, using a different synthetical argument. 
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6. — Calculation of the average number of emitted photons. Second part. 


Let us now derive an explicit expression for the matrix element appearing 
n (23). 

Inverting (5) we have: 
(29) j(e) =| ax exp [— itkey}j(e) 


So we obtain 


(30) as" |9(k) |a) = = dar exp [— i(ka)]<as™ | jw) |a?> . 


Here, to introduce more explicitly the known behaviour of the external 
system, we pass to the Schrédinger picture coinciding with the Heisenberg 
picture at t=— oo. Then using (10) and (11) we have 


(31) <os"*|j(k)|ar) =| ae exp [— i(ka)] <a" (ao) | j(x) |a (#0) = 


= —— dæ exp [— i(ka)] US: (to) Tir(co) a;(0) | F(x) |; (@o)> 


where the operator j(x) is constant with time. 
Now, using (9), we have 


(32) fare exp [— i(ka)] <a; (to) | f(x) | os (o)> = As: (k) exp [— è(E;— E; — ko)to] è 
where we have introduced the quantities 
FE A (le) = [are exp [— th] lita). 
Substituting (32) into (31) we obtain 
(34) <as* 19(k) [ar >= 7 DA (fan U?*;s (to) Utr (to) exp [— «E; — E; — ko) do] - 


Here we remember that we have assumed that the perturbation causing 
transitions between the states |«;> is different from zero only during the time 
interval between 0 and 7. 

Then (15) and (16) hold, and subdividing the integration over the time 
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in eq. (34) in three parts: from — co to 0, from T to +00, from 0 to 7, we 
obtain 
0 


(35) Co" | (ht) | ory = ily A jy (k) oa) ay exp [— (HL, — E,— ky)a| + 


\ 7 
(ce) 
(Ce) 


+ = A ,(k) or) fas exp = E, ZA E, = ky) ao | SE 
LA E 
L 


+ 2 A a | de, US,(c0) Us (to) exp [— (4, —E, — ky) = 


a) 


0 
asa (54, (k) Uk; (7) (250, — E, i o) + E, na E, mi 4 an 
Bi DA ,(k) U:,(T) exp [— E, — E,— ko)T] (now, — H,— k,) — E te 1 JE 
i 44 y Kg 


T 


on > Ay; a fa U;s(%0) Ui, (40) exp[_UE,—E,—k) à) : 


0 


Substituting this result in (23), a general expression for N,, is obtained, 
which may be calculated when we know the eigenstates of the unperturbed 
external system (energies Z, and matrix elements A,,(k), formula (32)) and 
the transformation matrix U‘(x,) of the interaction picture as a function of 
time. 

If the external system is constituted of electrons and positrons treated 
with second quantization, then we have: 


(36) (a) = 1p (x) yy(x) 
and (33) becomes: 
(37) A,(k) = iefars exp [— tk-x] <a; | p(x)yy(x)|a,>. 
Let us consider the case in which the external system is constituted of 
one electron in an external field. Calling |f,) its eigenstates in the configu- 


ration space, and w;(x')= <a'|f, its time independent unperturbed eigen- 
functions, we have in this case from (37) (1°): 


(388) A,,(k) = ie ara exp Fa ik:x/]y;(x')yw(x") = iecB;lya exp [— ik-x]y |B , 


where here x is an operator. 


(5) Or putting j(x)= y d(x —x’) into (32). 
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7. Proof of the finiteness of the average number of photons emitted in a 
transition between two bound states. 


Let us now discuss the conditions for the finiteness of N,,. 
Calling dQ, the element of solid angle in the k E the expression (23) 
may be written more explicitly 


(39) Fr pce SAL xf a, ko | Ca" | jo (Ie) [aa [2 

Now we remember that we are treating the case of emission of photons 
the energy of which is sufficiently low to allow this kind of approximate treat- 
ment. So the integrals in (23) or in (39) are extended only to a finite region 
V, of the k space. If this region does not contain the point k= 0, N,, is cer- 
tainly finite, otherwise an infinite amount of energy would be emitted: in fact 
an infinite number of photons of finite energy would be emitted. 

Therefore a divergence of the integrals appearing in (23) or (39) may only 
occur in the neighbourhood of k= 0. 

So W,, is finite when the integral over (k) in (39) does not diverge in the 
neighbourhood of k= 0. 

A sufficient condition to be in this case is that the matrix element 
<a" |4(k)|a'" be a steady finite function of k, at kj) = 0, for all the direc- 
tions of k. 

A more general sufficient condition is that 


(40) CR) loc? > = f(k) + g(k) O(ko) , 


with f(k) and g(k) steady finite functions of k, at k,— 0, for all the directions 
of k. In fact such an expression gives rise in (39) to integrals of the form 


È dkykyd(k) and i dkokod? (ko) 5 


which are finite (19). 

Now we shall prove that N,, is finite when the photons are emitted in a 
transition between two unperturbed energy eigenstates |x,) and |«,> be- 
longing to a discrete spectrum. We firstly suppose that no other unperturbed 
energy eigenstate with the same energy of the states |x,) or |«,> occurs. 


(16) For the meaning and the value of the second integral see footnote (°) of re- 
ference (4). 
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Then we may write again formula (35) in the form: 


(41) Cas" | i(k) |a) = = Us: (T)(A,,(k) — Ass (k) exp [th T]) + 

+ im 6(Kko) Us,(L) (A+, (k) + A,s(k) exp [iko T]) — 

e 1 e 3 
De (54,0 U;;(T) PIERRE a si (E) Us (7) 

i 1 Wa: 2 
exp i —B, —he)T] pp) +(Z A4 (MU: (ME, E, hy) + 
+ > A.:(k) Ui,(L) exp [— i(#; — E, — k)T]O(E—E.— to) + 
ifs 


+ | dx exp [— i(kar)] as" (20) | j() [22° (@0)> 5 


0 


where the last term has been written in its original form (31). Now we show 
that (41) has the form (40) in the neighbourhood of ky = 0. 

We observe firstly that, owing to the assumption cf the discrete spectrum, 
the fourth term of (41) is 0 in the neighbourhood of k, = 0. 

For the same reason the factors 1/(E,— H;—k,) and 1/(#;— E,—k,) in 
the third term of the right hand side of (41) are always finite in the neigh- 
bourhood of k,=0. Therefore each term of the sums over j and 7 is finite 
and it is easy to show that even the result is finite for the same values of ky. 
A rigorous proof is given in the Appendix II. 

The last term in (41) is certainly finite because the integration over the 
time is extended to a finite interval, and the integrand is finite owing to the 
supposition of bound states (17). 

At first sight the first term of (41) could give rise to an infinity of the 
form 1/k,. But this docs not happen because the last factor of this term is 0 
for tke— 0: 

In fact we have from (33): 


(42) 4,10) [da |j(2) |) - 


ra 


Now <«,|j(x)|a,> is the average value of the density of current in a statio- 


(©) A direct proof may be given starting from the explicit expression (35). Then 
the finiteness of this expression may be proved using formula (II.3) of Appendix II 
and arguments analogous to those used in Appendix II to prove the finiteness of the 
third term of (41). 
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nary state, so that (15), 


(43) div Kar | (x) | a,» = 0. 
Therefore a vector field I(x) exists, for which 


(44) <a, |j(x)|a) = rot I(x) . 


Then we have 


(45) A,, (0) — [as rot I(x) . 


The vector field I(x) may be chosen so that it tends to zero as |x| tends 
to infinity. In fact <x,|j(x)|x,) decreases exponentially as |x| tends to infi- 
nity, owing to the assumption of |x,) being a bound state, so that (44) be- 
comes rot I(x) ~ 0 in the neighbourhood of the infinity of |x|, and a solution 
of it is I(x) — 0 in the same region. 

Now let us consider a component of the vector A,,(0), for instance the 
third one A,,g(0). It is: 


OI,(x) 0I,(a) 
OX, 043 


(46) Ayr (0) = | dal | - | do (n, (x) I5 (a) — ng (x) L(x)) , 


o 


where we have transformed the integration over the volume into an integration 


(18) Indeed let us consider the unperturbed problem for the external system. The 
electric current in Heisenberg picture j%x) in this problem satisfies the continuity 
equation : 


__ Ojo (a) 
Cay © 


div j°(x) = 
So we have: 
5 E 0 : 
div Ko, | J° (x) |,» = Mor (CA (x) |) E 
Co 


Then, passing to the Schrôdinger picture, which coincides with the Heisenberg picture 
at #= — co, and using (9), we have 


5 d | d :; 
div (a, | j(a) |a,> = wr Co (00) [Fo (x) |, (20)) = pEr Car |jo(x)| a) =O. 


Now at t= — co the perturbed external system coincides with the unperturbed ex- 
ternal system, so that the result holds also for the perturbed external system and it 
may be written in the form (43). 


1116 


ON THE EMISSION OF LOW ENERGY PHOTONS FROM A QUANTIZED SYSTEM 205 


over a surface at the infinity, n(x) being the unit vector normal to the sur- 
face element. 

Then we see that the result is certainly zero because I(x) goes to zero 
exponentially as |x| tends to infinity. 

Therefore it is A,,(0)=0. 

So the first term of (41) has the form 0/0, as k, tends to zero, and it is easy 
to see that it tends to a finite value, as k, tends to zero (19). 

The factor which multiplies the 6(k,) function in the second term of (41) 
is certainly finite for k,=0 (we have shown that it is zero), so that (41) has 
the form (40). 

In this manner we have proved that the number of photons N,, emitted 
in a transition between two discrete unperturbed energy eigenstates |x,) and 
|a,> is finite in the case in which no other unperturbed energy eigenstate with 
the same energy as the states |a,> or |«,> exists. 

It is easy to extend the proof to the case in which there are other discrete 
unperturbed energy eigenstates with the same energy as the states |x,) or 
|a,>. This is done in Appendix III. 

The case in which the unperturbed energy of the external system has a 
continuum spectrum containing also the energy of the state |a,) or |«,> is 
more complicated. The extension of the proof to this case is treated in Ap- 
pendix IV (5°). 


8. — Average number of photons emitted in a transition between two states 
belonging to a continuum spectrum. 


It is easy to modify the formula (35) for the case in which the initial and 
final states of the external system belong to a continuum spectrum. Let p 
be a set of parameters labelling the unperturbed eigenstates |x,), of energies 
_#H,, which belong to the continuum spectrum. 

Let o(p) be the weight function of such a representation. Then, instead 
of (35), we have, for a transition between the states |x,) and |a) of the 


(19) In fact from De l’Hospital’s rule we have: 


Now this limit exists certainly and it is finite because 4,,(k) is the Fourier transform 
of a function of x which falls to zero exponentially as |x| tends to infinity. 
(29) The notation used in Appendix IV is explained in the next Section 8. 
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(47) oe |j(R)]ar) = 


= i [arena rik) Ua (D (mote, E, — ko) + (a i e) + SRE 
J SR cache melee 


+ | po) A, (#6) Up-(2) exp [—i(B, — B,—h) 
| Al 3 À bn At ; SAAS 1 + È 
it 6(#,— E,— ky) — ron e) + à st; +) dpap'o(p)o(p')Ay,(k) 
‘ 


; Jar Ur. (Lo) Ur: (10) exp [— (En — Er — ky) Xp] + > È ; 


ù 


where the terms > ...; >... >... have the same form as in (35), and the 
7 DI tj 

summations must be extended over the eventually existing discrete unper- 

turbed eigenstates. 

We see from this formula that the arguments used to prove the finiteness 
of N,, in the case of a transition between two bound states do no more hold. 
This fact is due to the terms containing the factors i/(H,— H, —k)) and 
1/(H, — Es— ko), which increase arbitrarily for kj = 0 when p has values in the 
neighbourhood of r and s respectively. 

Therefore the most important contributions to (47) come from these two 
terms for values of the variables of integration p near r and s respectively. 

Hence assuming that U{,(7) is continuous for p= r and U$,(7) is conti- 
nuous for p=s, we have, neglecting the terms of (47) which remain finite 
as ky tends to zero: 


a ae 1 1 
(48) co" lho) = — Sef ap ot p)( Ave SLE MAIN ETRS a) 


The approximations used to derive this simplified formula are the same 
used by Res Jos (®) and by JAUCH and ROHRLICH (?) in the calculation of 
W”. in the particular case of the transition between two plane waves of one 
charged particle. 


®. — Case of the transition of one electron between two definite momentum states. 
The expression (48) may be easily calculated in tke particular case of a 


transition between two plane wave states of one electron. The explicit cal- 
culation is given in Appendix V. 
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The same result may be derived also using the general property (28). 
Indeed the first order perturbative calculation of the probability amplitude 
of emission of one electron in this particular case has been made by RES 
Jost () and by JAUCH and ROHRLICH (?) with the same approximations; from 
the result it may be deduced: 


74 e? (pre) (pe) |? 
49 We ar og) | È 
ra 2 SAI pa) (0.5) 


Now, from (28) it follows that the same expression gives the average number 
of emitted photons N,,. This result is confirmed from the direct calculation 
of Appendix V. 

The expression (49) coincides with the expression found by THIRRING and 
TOUSCHEK (*) for the average number of photons emitted in the analogous 
classical problem, with analogous approximations. 

Therefore the number of photons emitted in the transition between two 
definite momentum states is the same in the classical and in the quantum- 
mechanical problems, provided the duration and the particularities of the 
transition are neglected, as it was made deriving (48) from (47) (23). 


10. — Extension to the case of a transition between two non-stationary states. 


I observe that the arguments leading to formula (18) for N,, hold also 
when |x,) and |x;) are non-stationary states of the unperturbed external 
system: in this case «is a complete set of commuting observables having |x,) 
and |x,) as eigenstates. Moreover, all the arguments leading to formula (23) 
hold also in this case, so that the average number of photons emitted in the 
transition between two non-stationary states |f,) and f,> is given by: 


a if 1 
50 ORE SERRA 
Da [SCE Cai 


| dk (Re?) | CB |.) | BD |? - 


Now let us express |f,) and |f,) as a superposition of stationary state. 
|x;). Then it is 


(51) che \j(e) | se) = > <Bs lea Cm [Bn Cani | i(k) | om, a 


(2) So this rigorous proof confirms the arguments given by JAUCH and ROHR$S 
LICH (2) to deduce W,, in this particular case. 
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So this problem too is reduced to the calculation of matrix elements of 7(k) 
between stationary states, and may be solved using formula (35). 

Moreover from this remark we deduce easily the finiteness of the number 
of photons emitted from an arbitrary system-in a transition between two states 
that are superposition of bound states. 


11. — Remarks on the case of a transition caused by a constant perturbation. 


Till now we have studied the problem of the emission of low energy photons 
in a transition between two states, caused by a perturbation of duration limited 
in time. 

There are in general conceptual difficulties to extend the problem to the 
case of a transition due to a time-independent perturbation. 

Let us first consider the case in which the effect of the time-independent 
perturbation on each particle is not limited to a finite time interval. 

Then the definition of a transition probability per unit time in the case 
of a time-independent perturbation implies the existence of times that are 
long compared to the periods involved in the transition, and small compared 
to the average life time of the initial state (2?). 

Such a condition is certainly not satisfied in the problem of emission of 
photons of sufficiently low energy in a transition between two discrete states, 
caused by a time-independent perturbation. 

Indeed for sufficiently low energy of the photons the periods of the photons, 
the emission of which is studied, become longer than the average life time of 
the initial state. 

So the problem of emission of photons of arbitrarily small energy in a 
transition between two states caused by a constant perturbation has in this 
case no meaning. 

The extension of these considerations to the case of a transition between 
states of the continuum, however, requires further examination. Indeed in 
this case the effect of the perturbation on each particle may have a finite du- 
ration, even when the perturbation is constant in time. 

Here we make only some remarks about this problem. 

Let us consider a transition between plane waves due to a time independent 
perturbation (usual S matrix problem). 

Then we have to distinguish whether the perturbation energy is different 
from zero only in a finite region of space or not. 


(22) See P. A. M. Dirac: The Principles of Quantum Mechanics (Cambridge, 1947), 
page 181. 
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In the first case (S matrix problem with external field) the effect of the 
perturbation causing the transition on each particle is limited to a finite time 
interval. It is then possible to define the problem of emission of photons of 
arbitrarily low energy in a transition due to the external field, and we can 
apply the treatment of Sections 10 and 11, which in particular leads to the 
Bloch-Nordsieck divergence. 

In the second case (S matrix problem without external field, for instance: 
relativistic quantum-field theory) the effet of the perturbation on each particle 
is not limited to a finite time interval. 

Therefore it seems reasonable to conclude that in this case the problem 
of emission of photons of arbitrarily low energy in a transition due to the 
interaction of the fields of the system has no meaning. Therefore the Bloch- 
Nordsieck divergence does not occur in an S matrix problem without external 
field, because the problem itself to which it refers does not exist in this case (?3). 


12. Remarks on the finiteness of N,, in real experiments. 


In Sections 7 and 10 I have proved the finiteness of the number of photons 
emitted from an arbitrary quantized system in an arbitrary transition between 
two bound states or between two superpositions of them. 

This result is sufficient to conclude that in all the real well determined 
experiments the number of photons is finite. In fact it is reasonable to assume 
that in all the real experiments the initial and final states are bound states 
or superpositions of them, otherwise an infinite region of space would be 
necessary to perform the experiment. 

This statement seems at first sight to contradict the fact that most of the 
experiments are undertaken with the purpose of measuring transition proba- 
bilities between plane wave states. However a further analysis of the experi- 
ments shows that only non-stationary states are prepared, that during some 
time interval behave as like as possible as plane waves. But before and after 
this time interval the particles behave in quite a different fashion, and generally 
sooner or later they come to rest in some part of the laboratory. Therefore 
the consideration of the experiments as scattering experiments is only an ap- 


(28) Therefore, if we start from the knowledge of the S-matrix without external 
field, the Bloch-Nordsieck divergence is an interference effect that occurs when a super- 
position of plane waves is considered, so as to form ingoing and outgoing wave packets. 
Indeed in this case the effect on each particle of the perturbation causing the transition 
has a duration limited to the time interval during which the wave packets interact 
‘on each other: so the consideration of wave packets instead of plane waves has an effect 
similar to the effect of an external field. 
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proximation, which is very good when the emission of photons of sufficiently 
small wave length is considered, but becomes inappropriate when we consider 
the emission of photons the wave length of which is of the same order as the 
dimensions of the laboratory or greater. In this latter case the appropriate 
approximation is to speak of a transition between two bound states (24). 

It is important to remark that the finiteness of the number of photons is 
no consequence of the fact that every detector of photons cannot reveal photons 
of energy lower than some limit, but it is a property of the spectrum of the 
photon emitted from a source of finite size in the region of long wave lengths. 

Precisely it is easy to see that the finiteness of the number of emitted 
photons is equivalent to the property that the energy emitted per unit of 
frequency-range tends to zero as the frequency tends to zero. 


13. — Connection with the infrared divergence of the perturbative calculation. 


Let us consider the case of a transition between two definite momentum 
states. 

Then we know that in the perturbative calculation of the probability of 
emission of one photon (bremsstrahlung) a divergence of this probability for 
photons of low energy is found. Such a divergence is called infrared divergence. 

Now in this same problem even the average number of emitted photons 
diverges. This was well known in the classical case (Bloch-Nordsieck diver- 
gence) and I have proved it rigorously in the quantum-mechanical case, in 
Section 9. 

Then the application of formula (28) to this problem immediately shows 
that these two divergences, although they refer to quite different physical 
quantities, are the consequence of the divergence of the same mathematical 
quantity, and therefore they are intimately connected to each other (?°). 


14. — Conclusions. 


The principal results of this work may be so summarized: 


1) I have given a general treatment of the problem of the emission of 
photons from a quantized system, when the reaction of the photons on the 
system is neglected. Till now such a problem had been solved only in the 
case of emission from a classical system. The most interesting quantity is the 


(24) See also (5). 
(25) See also (5). 
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average number of emitted photons, which is given from the general formula (23) 
or, more explicitly from (35). 


2) I have proved that the average number of photons emitted in an 
experiment, which can really be made and therefore need not an infinite region 
of space, is always finite: this is a particular property of the emission spectrum 
in the neighbourhood of the energy zero (Sections 7, 10, 12). 


3) I have proved rigorously that the emission of low energy photons in 
a transition between two definite momentum states of a quantum-mechanical 
particle is the same as in the case of a classical particle, provided the duration 
of the transition is neglected (Section 9). 


4) A connexion has been found between the results of the perturbative 
calculation and those of the exact calculation, which holds exactly when the 
reaction of the photons on the charges is neglected; the average number of 
photons emitted in a transition is equal to the first term of the expansion in 
powers of e of the emission probability of one photon in the same transition 
(Section 5). 


5) The application of this property to the case of a transition between 
two definite momentum states leads to the conclusion that the Bloch-Nord- 
sieck divergence and the infrared divergence of the perturbation calculation 
have their origin in the divergence of the same expression, notwithstanding 
they are divergences of different physical quantities (Section 13). 


APPENDIX I 


Let us call W,,, the probability of emission of N photons in a given mo- 
mentum range, caused by a transition of the external system from the state |«,»> 
to the state |x,). Let N,, be the average number of photons emitted in the 
same momentum range and in the same transition. 

Then it may be asserted without need of direct calculation that in our 
approximation W,,, is given by the Poisson distribution: 


—y 


| N Ts x 
(1.1) War wi XP [N]. 


In fact let us consider all the processes in which a photon is emitted in 
the infinitesimal energy range dk,, while eventually other photons are emitted 
in other energy ranges. Then we observe that each photon cannot interact 
with the other photons in this approximation, because we neglect the reaction of 
the photons on the charges, and photons can interact on each other only through 
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charges. So the emission of a photon in the energy range dk, is not affected 
from the emission of other photons in other energy ranges in the same tran- 
sition. Therefore the emission of a photon in a given infinitesimal energy 
range is an event independent on the emission of other photons, so that we 
may define an emission probability of a photon per unit energy range, and 
one deduces in the usual way that the probability Wy,, of the emission of 
N photons is given by the Poisson distribution (I.1). 

So the knowledge of N,; is sufficient to individuate any transition probability. 


APPENDIX II 
We show here that (see Sect. 7): 
(HET) D (A; (k)e) UT) > 
î rg o HD; 
is finite in the neighbourhood of k, = 0. 
From (33) we have: 


(12) XI, = dl (| fato exp ik-alj (x)]a) |? = 


= CO, | 


( fare exp[— ik: ‘a1j(2))-el 


So the sum Y|(4 jr(O)e) |? is equal to the expectation value in the state |a, 


of the square of a component of the total electric current |d*xj(x). 

Now, if the theory has a meaning, the expectation value of this quantity 
in a bound state must be finite (**). So the expression (II.2) must be finite 
for k — 0, and it is also finite in the neighbourhood of k — 0 (?). 


(26) The fulfilment of this condition can be verified directly when the actual ex- 
pression of j(x) is given. 

For instance in the case of one Dirac electron we have from (38) (x is here an 
operator) : 


>| (Air(k)e) [P= 2118 [—tk-x](ve)|B,>|? = <Br](ve)(ve)]B,) = <Br|Pr> = 1: 


(27) Indeed it may also be written: 


par ry t= | aba exp [— ik-(a — x’) Kay [fe (æ)je (%') op) - 


So it is a Fourier transform of a function of x which falls to zero exponentially as 
|x| tends to infinity, owing to the assumption that |«,> is a bound state. 
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So 
(11.3) 2! (A slk)e)|?< Zi (A,(k)e)|}=A?, 
where A is finite. 

Moreover it is 
aS) SS DT. 


si 


Let us call AE,» the smallest of the quantities | E, — E;|, j<r. Let us 
suppose ky << AÉ,m: Then from (II.4) it follows 


| 1 2 1 
(II 5) > U (7) E, = E; — ko (A rim ad E) 


From (II.3) and (II.5) we have, using the Schwarz inequality: 


È He A 
ead a 4 sr (K)e) Us(1) E,— E;— ko n AEG, Sete ‘ 


aAr 


Therefore the expression (II.1) is finite when kj << AZ: 


APPENDIX III 


Let us suppose that the initial and final states |x,,) and |x,,) belong to 
multiplets |«,,>, |æ,), 4, j=0, 1, 2, ... of bound states, where the vectors in 
each multiplet have been chosen so as to be orthogonal each to other. 

Then the two first terms of (41) become 


(111.1) ae UL) A, (E) > Ut, (T) Aus, (k) exp 0) Ti 


+ i d(Ko) (EX Vig, (L) Arg (k) +S Usro(P)As,(k) exp[ik,T]) , 


whereas the remaining terms have the same form as in the case of non- 


degenerate states. 
The argument given in Sect. 9. proves even in this case that A,,,,(0)=0 


and A,,,,(0) = 0. 
However it is also 4,,,,(0)—0 and 4,,;(0)=0. 


14 - Il Nuovo Cimento. 
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In fact from (33) we have: 


(111.2) Ar (0) = | "Car, |j(#)|xp 


Now 
ŒIL) <0, f(x) a, = 
= (Ka, + op, | G(X) le, + Or,» = <a, |j(£)|x,) — on |F(%) ar )) - 


The only properties of |«,> used in Sect. 7 to prove that [dax lix) |= 
are that |x,) is stationary and bound. Now also <«,,+-«,,| has these two pro- 
perties, therefore : | 


fate can + he |XX) |a + an) = 0, 


and from (III.3) it follows 4,,, — 0. 
So all the arguments given in Sect. 7 may be used to show that also in 
this case the expression (41) has the form (40) so that N,,, is finite. 


APPENDIX IV 


Let us suppose that the external system possesses a continuum spectrum 
having energy eigenstates equal to the energies £,, E, of the initial and final 
bound states. Then in expression (41) the first two terms and the last term 
conserve their form, whereas the third and the fourth terms must be written 
with an integral besides the sum; so the following integrals must be added 
to the expression (41): 


7 1 
(IV.1) [remi MM lu 
e 7 1 ii 
+ | dp o(p)A s»(k) Uz, (T) exp [— «E, — i, he) pred ko na 


Br i dp 0(p) Ar (k)U!:(D)0(E;— E,— ko) + 
ae | dp A, (4) Uge(T) exp [— i(E,— E,— ky) T16(Ep — E,— ky) - 


The symbols used here are defined in Sect. 8. 
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Now let us choose the energy £ as one of the parameters p labelling the 
continuum eigenstates, and let q be the remaining parameters. 

We suppose that the parameters g are chosen so that o(q, E) is a con- 
tinuous function of HL. 

Then the first integral of (IV.1) becomes 


1 
(IV.2) foafazon, EF) Ag, mr(k) Usq.2(L) i, . 


This expression is certainly finite, for the argument of Appendix I, when 
an arbitrarily small domain H,—k,—e< H< E,— k, + 8 is omitted in the 
integration over E. It remains then to prove the finiteness of 


Er-ko+€ 


n aft 
(IV.3) IL À dE o(E, q)Aw.rr(k) Uan) FER 


E,-ko-€ 


where e is arbitrarily small (but finite). 
Now (IV.3) is approximately equal to (?§): 


Er-ko+€ 


(IVA) | ina | dr 


Er -ko-8 


E: E ki 


Here the last integral must be interpreted as a principal value because 
the factor 1/(E, E — k,) has been obtained from the formula 


a 


DO (Eko = RE, — E-k)+7—5_x' 
corn — Ko 


so that the expression (IV.4) is approximately zero, and therefore it is cer- 
tainly finite. So the expression (IV.3) too is certainly finite. 
The third integral in (IV.1) may be transformed analogously into 


(IV.5) Jaufar 0(4, BE) Arr (k) Usa.» (T)0 (B,— E — ko) = 
= dgo(g, E,— ko) Aa 2,-xpr (K) Usa.z;-t9 (7) ; 


and it is certainly finite for the argument given in Appendix I. 


(28) This statement uses the facts that Us, (7) is certainly a continuous function 
of E and that it is very reasonable to assume that also 4,,,,,(k) is a continuous function 
of E. 
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The application of the same arguments to the second and fourth integrals 
of (IV.1) is trivial; so the whole expression (IV.1) is finite. 


Therefore even in this case the number of emitted photons - “Vena to 
be Mara: uy 


APPENDIX V 


Let us apply the expression (48) to the case of a transition between two 
plane waves of momenta p,, Ps. 


Let us use the momentum p to label the unperturbed eigenstates. 


Let the wave functions be normalized to one particle per unit volume, 
so that o(p) =1: 


Wp(x) = u, exp [ip-x]. 


Then (38) gives: 
ae ie) iel ty) | a2 exp [i(p:— p;— k)-x] = te(up,yup,)0(pi— pi— E). 
Now (Up,yUp,) is a continuous function of p;, p;; so for small ky it is 


Ap,p;(k) = 1¢(Up,7Up,)0(p à — ps — k) = ef 0(pi— pi— k). 
Therefore (48) gives 


(V.1) ato | i(k) ) | oe; DI — GU): 


Pr fas Ds (qs ee -| 
(E fe pÒ(Py—P Pier ee Al pÒ(p — ps ) 7 ee = 
ge (Dr al Ps Jr 
Wi DCR ee er Erra 


Now for small k, it is: 


-k 
LES = ah he 


Pr 1 Ik 
E,— HSE ko = Ee A ko = È (p,k ss Eko) = (pk) 
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so that (V.1) becomes 


<a | j(Ie) 


ds = es, | Pr ook 


Substitution into (23) gives (49). 


RIASSUNTO (*) 


Risolviamo qui il problema dell’emissione di fotoni da un sistema generale quantiz- 
zato quando si trascuri la reazione dei fotoni sul sistema. Finora tale problema era stato 
risolto solo per il caso dell’emissione da un sistema classico. Deduciamo un’espressione 
generale per il numero medio di fotoni emessi in presenza di una transizione tra due 
stati arbitrari del sistema prodotto da una perturbazione dipendente dal tempo. La 
discussione di questa espressione porta a concludere che il numero medio di fotoni di 
bassa energia emessi in qualsiasi esperimento realmente eseguibile è finito: questa 
conclusione generalizza al caso di una sorgente quantizzata un risultato precedente- 
mente ottenuto dall’autore per il caso di una sorgente classica. Inoltre, trascurando 
la reazione dei fotoni sulle cariche, dimostriamo qui che il numero medio di fotoni 
emessi è uguale al primo termine dello sviluppo in serie di potenze di e della proba- 
bilità di emissione di un fotone. L'applicazione di questa conclusione al caso della 
transizione tra due stati d’impulso definito mostra che la divergenza di Bloch-Nordsieck 
e la divergenza infrarossa del calcolo perturbativo hanno la loro origine nella divergenza 
della stessa espressione matematica, per quanto abbiano un significato fisico completa- 
mente differente. Finalmente dimostriamo che il problema studiato qui da un punto 
di vista generale non ha significato nel caso di una perturbazione indipendente dal 
tempo. 


(*) Traduzione a cura della Redazione. 
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The Worldwide Distribution of the F2 Layer Electron Density: 
Seasonal and Non-Seasonal Variations 
and Correlations with Solar Activity. 


F. MARIANI 


Istituto Nazionale di Geofisica - Roma 
Istituto di Fisica dell Universita - Roma 


(ricevuto il 20 Dicembre 1958) 


Summary. — The worldwide behaviour of the twelve-month period 
(seasonal and non-seasonal) variation of the maximum electron density N 
in the F2 layer appears noticeably different at noon and at midnight. 
In each case, it is controlled by the geomagnetic field and by the solar 
activity but, whilst at noon it shows large asymmetric features in the 
two hemispheres, at midnight, on the other hand, it is fairly symmet- 
rical. It is not possible to account for the experimental results by assuming 
some external source of radiation producing non-seasonal effects. One 
is led to assume some terrestrial cause of asymmetry: for example, some 
effect of general circulation in the upper atmosphere. With respect to 
the correlation of electron density with solar activity, there appears 
again some noticeable asymmetry in the two hemispheres; in particular, 
the northern hemisphere appears more influenced by the solar hydrogen 
filaments, especially for increasing latitude. The explanation of these 
asymmetries, in terms of some direct influence of the sun, is also puzzling, 
so that, in this case, one is led to think of some more or less unknown 
effect of solar corpuscular (or electromagnetic) radiation on the upper 
atmosphere, for example on its movements or its conductivity. At the 
present stage of our knowledge, we cannot exclude the possibility of a 
similar cause for both the asymmetries of the twelve-month period vari- 
ation and of the correlations with solar activity. Further observational 
evidence is required. 


1. — Introduction. 


A very important question in ionospheric physics is the study of the pe- 
riodic and secular variations of the electron density, particularly those of the 
2 layer which exhibit many interesting anomalies. Many authors have con- 
sidered various aspects of these anomalies, but the problem of the seasonal 
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and no#-seasonal variations of the electron density in the F2 layer, first con- 
sidered by BERKNER-WELLS (1) and ECKERSLEY (?), requires a worldwide sta- 
tistical investigation, to elucidate their dependence on geographical or geo- 
magnetic latitudes, on solar activity, ete. 

Such a statistical investigation is now possible, since we possess a nearly 
continuous collection of experimental data for many ionospheric observatories, 
adequately distributed over the Earth, and for time intervals sufficiently long 
to warrant a more definite experimental basis for future theoretical investi- 
gations. 


2. — Experimental data and method of analysis. 


We have used the ionospheric monthly median values of foF2 at noon and 
at midnight, taken from « Ionospheric data » volumes of the National Bureau 
of Standards; in our opinion, the use of median values is preferable to that 
of the mean values, because they are less influenced by perturbations; and 
we are principally interested in the « normal » behaviour of the F2 layer. With 
regard to the solar data, we have used the Wolf number À and sunspot areas 4,, 
respectively measured in the observatories of Zurich and Washington; for 
the other solar parameters we have used the values measured in the Astro- 
physical Observatory of Arcetri. 

For the analysis of the data we have used Vercelli’s selective method of 
periodogram analysis (VERCELLI (*4)), which is based on symmetrical linear 
combinations, so that it does not alter the effective phase of every periodic 
component, while the corresponding amplitude is reduced for a known amount, 
depending on the coefficients of the linear combination. Thus one can obtain 
the correct phases and amplitudes. We have first applied the periodogram 
analysis to the series of median values of N= (foF2)?, having first obtained 
the smoothed time series and, afterwards, the series representing the periodic 
component N,,, having a principal period of twelve months, and the secular 
component N, having the eleven-year period of the solar activity. 

In effect, in the scheme we have used, together with the principal period 
of twelve months, there remain in the series of N,,, with amplitude compa- 
rable with that of the principal component, also the periodic waves of eight 
and ten months, which however, except the cases of Delhi at noon and of 
Singapore at midnight, will result practically absent. With the same method 
we have obtained the secular series for the solar parameters. 


L. V. BERKNER and H. W. WeLLS: Terr. Magn., 43, 15 (1938). 
T. L. Ecxerstey: Terr. Magn., 45, 25 (1940). 
F. VERCELLI: Ric. Scient., 51, 364 (1934). 

F. VerceLLI: Memoria ccLxxxv of the « R. Comitato Talassografico » (1940). 
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We note that the values of N at noon may be considered approximately 
proportional to the ionization intensity just at the level of maximum electron 
density of the F2 layer. 

Data from about thirty observatories have beens used, as indicated in 
Table I, which gives their geographical and geomagnetic coordinates and their 


MB I (6). 


Geomagn. coordinates Geogr. coordinates Macnetic 
Observatory di x 
Lat. Long. Lat. Long. I 
Resolute Bay | 82 N 289 E 75 N 265 E 85 N 
Reykjavik 70 | 71 64 338 76 
Point Barrow | 68 241 71 203 80 
Kiruna ler en 116 68 20 = 
Anchorage 61 258 61 210 73 
Lindau 52 94 52 | 10 66 
Freiburg | 50 90 48 | 8 64 
Washington 50 350 39 283 71 
San Francisco | 44 298 37 238 62 
White Sands 41 316 33 259, 62 
Baton Rouge 41 334 30 269 63 
Wakkanai 35 206 45 | 142 60 
Tokyo 25 206 36 140 50 
Yamagawa 21 198 sl 131 45 
Maui 21 | 268 21 203 38 
Delhi 19 149 29 at 43 
Bombay 10 144 19 73 26 
Madras 3 150 13 80 lt 
Guam 3 212 13 145 14 
Huancayo 0.6 S 354 E 128 | 285 E 4 N 
Leopoldville 3 84 4 15 33 S 
Singapore 10 173 1N 104 16 
Rarotonga 21 274 218 200 38 
Johannesburg 27 91 26 28 62 
Brisbane 36 227 27 158 57. 
Watheroo 42 186 30 116 64 
Canberra 44 225 35 149 65 
Christchurch 48 253 44 173 68 
Hobart 52 225 43 147 72 
Decepcion 52 7 63 299 56 
Port Lockroy 53 4 65 297 58 
Campbell Is. 57 253 53 169 76 
Terre Adelie 75 231 67 140 88 
(*) The data of this table, sligthly different than those given by MARIANI (*) are desumed 
from the data collected by the C.S.A.G.I., now at our disposal. Moreover the table gives the 
correct coordinates of Christchurch which there, for an error of transcription, are mistaken. 
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magnetic inclinations. The experimental data at our disposal extend from 
January 1946 to December 1955. However the series of N, N,,, etc., extends 
from January 1947 to December 1954; because in the method of linear com- 
binations one loses part of the values of the starting series. Within this last 
interval we have selected, for the calculation of correlations, the interval 
May 1947—March 1954 which can be assumed as the complete decreasing stage 
of the last solar cycle. For some of the observatories indicated in the Table we 
have had at our disposal only shorter series of experimental data, so that the 
series of N,, and N are also shorter. 


3. — Seasonal and non-seasonal variations. 


31. Behaviour of N,,. — We have considered the experimental series of N 
at noon and at midnight: Fig. 1 shows the graphs of N,, at noon for the obser- 


N. 1948 1949 | 1950 1951 1952 1953 1954 


12 


/ \ a PTX Frans 
ASIN VA ANIA TI 
42017 UN VITA \y NR Sands (410) 


500 IN 


WA 


Wakkanai (370) 
Tokyo (350) 
Yamagawa (330) 
Maui (310) 


Delhi(290) 
Bombay (265) 
Madras (235) 
Huancayo (215) 


Rarotonga(185) 


Johannesburg (145) 


Brisbane (125) 


Watheroo (105) 
a (90) 


Christchurch (55) 
Hobart (25) 


Fig. 1. — In this Figure and in the successive Figs. 2, 3, 4, 5, 6 the numbers within 
parentheses indicate the ordinate corresponding to the zero ordinate for every graph. 
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ratories whose series of data are about complete for the decreasing stage of 
the solar cycle, arranged according to the geomagnetic latitudes. We note 
the following principal features: 


-1) the amplitude of N,,, for-symmetrical northern and southern lati- 
tudes, is clearly greater in the northern hemisphere; 


2) in the northern hemisphere, the times of maximum or minimum 
amplitude show some slight tendency to. be anticipated during the epochs 
of lower solar activity and at higher latitudes; 


3) in the southern hemisphere, as far as the latitude of Brisbane, the 
component N,, exhibits a substantial phase-agreement with respect to the 
northern hemisphere; however, during the years 1947-1948 of higher solar 
activity, the behaviour of N,, at Johannesburg and Brisbane appears as a 
transition stage to the behaviour at observatories at greater southern lati- 
tudes, which clearly show a phase-shift of 180° with respect to the other obser- 
vatories; on the other hand, with regard to the years 1949-1953, the times 
of maximum or minimum amplitude of N,, show a systematic tendeney to 
be anticipated during the epochs of lower solar activity and at higher lati- 
tudes; 


4) as a whole, the component N,, decreases for decreasing solar acti- 
vity; between the geographical latitudes of S. Francisco and Rarotonga, an 
appreciable temporary diminution of amplitude is noted around June 1948 
and June 1950; conversely, an increase of amplitude appears around June 1948 
at Watheroo, Canberra, Christehurch and Hobart. 


Concerning the behaviour of N, at midnight we refer to Fig. 2, from which 
we can infer the following principal features: 
1’) the amplitudes of N,,, for symmetrical northern and southern lati- 


tudes, are substantially the same; 


2') the times of maximum or minimum amplitude practically coincide 
with the solstitial months; 


3’) there is a substantial phase-opposition for symmetrical northern and 
southern latitudes; 


4’) as a whole, the component N,, decreases for decreasing solar activity. 


Some further interesting features, that one can draw from Figs. 1-2, are 
the following: 


5) the amplitudes of N,, are greater at noon than at midnight, in both 
the northern and southern hemisphere, with the exception, after 1949-1950, 
of the observatories at latitude higher than Brisbane; 
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6) the geomagnetic control of N,, either at noon or at midnight is clearly 
manifest; moreover there appears some «regional» behaviour, for example 
the noticeable attenuation of Ns at midnight in the observatories of Washing- 
ton, San Francisco and White Sands, which are «regularly » distributed only 
with respect to the geomagnetic latitude. 


Lindau (340) 
Fribourg(290) 
Washington (310) 
S. Francisco (305) 
White Sands (280) 
|Baton Rouge (260) 
A | Wak kanal (245) 
Tok yo (230) 


Yamagawa (215) 


Maui (190) 


Delhi (150) 


Huancayo (140) 


Rarotonga(110) 


Johannesburg(100) 
Brisbane (80) 


Watheroo (50) 


LR Li 
NV am ( D Canberra (40) 
NYE RY? eons 
AVAMD-AMY/ 


In order to obtain further information on N, for very low and high lati- 
tudes, we now refer to Fig. 3, where the graphs are arranged according to the 
geographical latitude of the observatories, or else directly to the monthly 
median values of fo#2. The graphs of Fig. 3, which refer to shorter time inter- 
vals, show the following features: 


7) at noon, for equatorial and high southern latitudes, the component N,» 
reach maxima and minima respectively at approximately the December solstice 
and the June solstice, with the surprising exception of Singapore, where the 
reverse is the case. 

For high northern latitudes the phase angle of N, varies almost continually 
with increasing latitude so that, at Resolute Bay, N» shows a «normal « be- 
haviour, with summer maxima and winter minima, thus concluding the prog- 
ressive phase shift already noted in the graph of Fig. 1. One notices a geo- 
graphical rather than a geomagnetic control of N, at high latitudes. 
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As regards the amplitude of N,,, it appears approximately the same, at 
corresponding equal southern and northern latitudes. 


1401 si Resolute Bay (140 
1 > Point Barrow (130) 
120 3 Kiruna (120) 
| Reykjavik (110) 
100 i Guam (85) 
Anchorage (100) 
80 4 | 
Ke BAS ge Singapore (55) 
A | AVES Leopoldville (25) 
20 | 
| Port Lockroy (0) 
(0) —— 2} 
160 + ZA i Bay (160) 
IN ge alee E Kiruna (150) 
140 = | 
120 | | Oy eM ewe UN ua (130) 
G (100) 
100 uam (1 
80 Singapore (75) 
Leopo/dville (60) 
ii rei 
60 Midnight | 
| Decepcion (45) 
e 7 Port Lockroy (30) 
20 
0 
1950 | 1951 1952 


Fig. 3. 


8) At midnight, the N,, shows phase-agreement at high and equatorial 
southern latitudes, again with the exception of Singapore, where a clear twelve 
month variation is not apparent. In respect to high northern latitudes, Nj. 
shows «normal » summer maxima and winter minima. The amplitudes of N, 
are much larger at high southern latitudes than at high northern latitudes. 
The graph of Guam is very interesting: although both the geographical and 
the geomagnetic latitudes are northern, N,, shows a substantial phase-agree- 
ment with the values relating to southern latitudes; thus, if one assumes the 
line of phase-inversion of N,, as «ionospheric equator» this appears rather 
different from the geomagnetic equator (as well as from the geographical 
equator). Our definition of «ionospheric equator», which can be considered 
rather similar to that given by APPLETON (5), seems however to lead to rather 


(5) E. V. APPLETON: Journ. Atmosph. Terr. Phys., 1, 106 (1950). 
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different conclusions. Obviously a larger number of experimental data are 
necessary in order to derive the position of the ionospheric equator (or, perhaps, 
equators). 


Concerning the behaviour of the electron density for higher southern lati- 
tudes we have had only shorter time intervals of data at our disposal: how- 
ever, the visual inspection of monthly median values of foF2 for the obser- 
vatories of Terre Adelie (February-December 1951) and Campbell Islands 
(April 1950-October 1951) seems to suggest for N,, characteristics similar to 
those of Decepcion and Port Lockroy. Another interesting fact is the syste- 
matic occurrence, for foF2, of values rather higher at midnight than at noon, 
during the winter at the observatories of Decepcion, Port Lockroy; Terre Adelie 
and Campbell Islands; no indication of any similar fact appears for the higher 
northern latitudes. 


3°2. Discussion. — As a first conclusion, the comparison between the N, 
relative to the various observatories, seems to suggest some interesting asym- 
metric features between the northern and the southern hemispheres, which, 
on the other hand, seem principally to concern the daylight side of the Earth. 

In order to attempt some explanation of the observed anomalies, one has 
to remember the results and the discussion of BERKNER and WELLS: they 
studied the variations of the maximum electron density in the F2 layer at 
Washington and Watheroo during a period of about three years and inferred 
the presence of a non-seasonal variation of the electron density, having a prin- 
cipal period indistinguishable from one year and an amplitude about equal 
to that of the seasonal variation; the amplitude appeared to maintain an 
approximately fixed ratio to the average background ion density. 

The authors showed that the non-seasonal effect could not be explained 
by the difference of latitude of the two stations from the equator, nor by the 
ellipticity of the Earth’s orbit; on the other hand, also the non-homogeneity 
of data could not account for the non-seasonal variation. As regards the source 
of the effect the authors were unable to make any definite suggestion; they 
did not exclude that, for example, the effect might be associated with sidereal 
causes, or with the varying position of the Earth’s magnetic moment with 
respect to the Sun, or with the circulation in the high atmosphere under the 
influence of the Earth’s magnetic field. Some other observational evidence 
on the non-seasonal variation is due to ECKERSELY, who suggested some even- 
tual association of the effect with JANSKY noise (JANSKY (°)). 


(5) K. G. JANSKY: P.I.R.E.. 21, 1387 (1933). 
(*) K. G. JansKy:]/P.I.R.E., 28, 1158 (1935). 
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At this point, our extensive results allow us to draw some definite con- 
clusions; for sake of simplicity we shall indicate as N°, and N°,, respectively, 
the non-seasonal (annual) and the seasonal components of N,,, the former in 
phase-agreement between the northern and southern hemispheres, the latter 
with a phase-shift of 180°. As is well known, it is still difficult to conclude 
whether the «normal» seasonal behaviour of the 7, electron density is the 
one with summer maxima and winter minima, or viceversa: for this reason, 
in the theoretical analysis of the worldwide variation of N,., we have con- 
sidered as «normal» seasonal variation the cases of both summer or winter 
maxima. 

At first, we have assumed a sidereal source or, at least, some association 
with some sidereal phenomenon to account for the non-seasonal effect. If it 
is so, we should have to expect that the effect depends upon the sidereal time: 
thus, if it is maximum at midday in December and January (i.e., midwinter, 
northern hemisphere), it should also be a maximum at midnight in June and 
July (iî.e., midsummer, northern hemisphere). In effect, if the behaviour of 
N,, at noon can suggest some sidereal dependence of the non-seasonal varia- 
tion, the behaviour at midnight clearly disagrees with such a hypothesis, be- 
cause in the southern hemisphere the component N,, is a maximum in De- 
cember and January, 1.e. in the austral summer; on the other hand, also 
during daylight, at very high northern latitudes, N,, reaches its maximum in 
June and July, as it does during the night. 

Moreover, any scheme of interpretation must give account of the phase 
opposition that N,, shows at latitudes higher than that of Brisbane (in parti- 
cular that of Watheroo, whose data were considered by BERKNER and WELLS) 
with respect to the northern latitudes, during the period of maximum solar 
activity; moreover, another singular aspect of the behaviour of N, is the 
gradual phase-shift depending on the latitude or on the solar activity, which 
is not clearly understood in the scheme of seasonal and non-seasonal varia- 
tions of N, both exhibiting solstitial maxima and minima. 

A last a priori objection that one can raise is the intensity required for 
the possible sidereal source of radiation, which would have to be of the same 
order as the solar radiation. However, one can assume the effective presence 
of two N°, and N°, components of about the same order, both reaching their 
maxima and minima in the solstitial months, and study their superposition 
in all the possible cases: we do not wish to labour this matter because, as can 
easily be verified, in every case, some systematic phase-disagreement or ampli- 
tude difference is obtained between the experimental and the expected be- 
haviour of N,, in the two hemispheres. 

As a first result of our discussion, the possibility of some important sidereal 
effect on the ionization of the F2 layer may be excluded. We can thus con- 
centrate our attention on solar or terrestrial causes: as we have seen, the non- 
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seasonal irregularities occur essentially during daytime so that in order to 
explain them we are led to study in some detail the relations between solar 
radiation and our atmosphere. 

In effect there is a non-seasonal variation of solar radiation due to the 
variation of the distance of the Earth from the Sun, which is a minimum at 
the December solstice and maximum at the June solstice: however, the dif- 
ference of solar intensity in the two periods does not exceed a few per cent, 
so that one is led to infer the effective existence of some asymmetric pro- 
perty of the F2 layer in the two hemispheres, principally or exclusively during 
daytime. We shall consider the question again in Section 5. 


4. — Correlation between electron density in the #2 layer and solar activity. 


41. First case: at noon. — We begin by considering the secular behaviour 
of N and its correlations with the Sun. For the parameters characterizing solar 
activity we have selected, together with the Wolf number #& and the sunspot 
area A,, also the areas of the hydrogen flocculi and filaments, respectively 
indicated with A, and A,, which are quiescent phenomena reappearing on 
the solar dise also for 
many consecutive ro- 
tations. 

If we plot the se- 
cular variations of N 
dior Ag, Ap 
(which we indicate - 


with a bar: N, R, etc.) 4 L Delhit130) 
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Fig. 4. — The areas 4, 

and A4y are the projected 

areas, expressed in 10-4 

of the solar disc; the 

areas A, are expressed 

in 1/(16.09)-106 of the 
solar disc. 
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5. — In this Figure and in the successive 
. 6,7 the arrows indicate the date: for ex- 
oi 1-50 indicates January 1950. 


The black 


points are the values of N for odd months. 


sponding relative maxima and an 
intermediate minimum in the nor- 
thern hemisphere; in the southern 


-hemisphere, instead, such « simi- 


larity » is not present with the 
exception of Huancayo and Ra- 
rotonga. 

Some better information is 
given by Fig. 5 and Fig. 6 which 
supply respectively the graphs | of 
N (R), and of 4,(R) and 4,(R); 
we do not give the graph of 
A,(R) which, except for some 
very small fluctuation, is well re- 
‘presented by the linear relation 
A,=16.09R, so that, as a con- 
sequence, one can infer that the 
behaviour of N(R) is identical 
with that of V(A,); a remarkable 
correlation between N and A, 
appears again and more clearly; 
conversely, no striking correlation 
between N and A,, appears, so 
that we shall disregard every 
dependence between N and Aj. 

Moreover, one can infer that 
the correlation between N and 


A, seems to increase for increas- 


ing northern latitudes; however, some increased dependence of N upon A, 


appears also at the observatories of Christchurch and Canberra. 


both in the northern and south- 
ern hemisphere, a sudden va- 
riation of the slope around 
the beginning of 1951 is appa- 
rent at about the same time 
as a similar sudden variation 
of the slopes of A, and 4,. 

In order to study the fore- 
going results, it may be obser- 
ved as a preliminary question 
that, in effect, the monthly 
median and mean values N 


Finally, 
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and À are referred to the conventional terrestrial month, so that they would 
have to be replaced by the analogous values N* and R* referred to the solar 
month; a calculation has, therefore, been made of the new series of solar month 
“median and mean values N* and 
R*, for Washington: the new 
graphs of N*(R*), as can be seen 
in Fig. 7, reproduce the features 
of the graphs shown in Fig. 5, so 
that one can infer that the featu- 
res of the secular variation of N 
have some substantial physical 
meaning. 

We have therefore assumed a 
multiple linear regression of N on 
R and A,, expressed in the form: 


(1) N=VN,1+6R+ y4,) 


and the simple linear regression of N on R and N on 4, 
(2) N=N1+«R), 
(3) N = N(1+04,). 


Table II shows the values of N,, 0, B, y, à, the simple correlation coefficients 
Tyr Tya, AN the partial correlation coefficients r,, Ap) THAgB (which will hence- 
forth be indicated by the simpler symbols #3; #13} 712.3; 713. attributing the 
indices 1, 2, 3 respectively to V, R, A,), calculated according to the standard 
linear correlation methods (OSTLE (°)); we also give, the standard deviations 
of the partial regression coefficients; finally, we supply the mean values %, 
Por Yo: do Of «, B, y, 6 for the eight observatories of the southern hemisphere 
and the corresponding standard deviations o,, 6,, o,, 0, together with the 
mean values Ao, Bo, 1%, A, for eight observatories located at symmetrical 
northern geomagnetic latitudes, which are indicated in italics. 

One sees in the Table that the coefficients «, 6, f have a statistical meaning, 
while this is not the case for all the coefficients y which cannot at first sight 
be considered significantly different from zero in the southern hemisphere 
(other than at Maui and Baton Rouge in the northern hemisphere): there 
might be some tendency of y to increase in value for higher southern lati- 


(8) B. OstLE: Statistics in Research (Ames, Iowa, 1954). 


15 - Il Nuovo Cimento. 


1141 


= : 9800 =*o | ogo =% |<? 2200 =o | |6s00=%2 | ; cui 


| 

i | gF00 =" 1610 = °¢ | et | gio =" | | 
| (961 | 606 | 296 | 666 | 200° L200" FI00 + 9FTO" E" | 100° + 970" FSE 2000070910" USE | Q1eqoH 
QT | GER | 096° | NS | FF00 + F800" | zoo cogo gst | 00 FS6F0" | LEE | SOOO FOFEO | EST | PMESE) 
Lacs | rés: | 906 | 166° | re00 Freoo | 2100 Frogo | see | 9100 680" | 268 | 000 FLFEO | OSE | BLEUE) | 
‘ | LEE | FEE | FOG' | CEG | LZ00°TEFOO | FI00 FCIS0 | LE | 9100 F FLO" | LSE | FOOD TSETO ELE | come | 
LOF | L96 | 896 | 266° | LI00 FSFOO | 6000 TESTO. LTE FI00 F0" ELE | 000 LOTO LEE eueqsug | 
IST OLO €96 | 866 | SI00 1007 9000 FSIO OSE  £I00 79080" LFF TOOOTIGIO CSE Smqsouuego f | 
get” | IR6 | 896 | L66° | TION FELON | 1000 79610 | 2°09 | LI00 TEFEO” [19 &000 T&FIO 9°09 eFnoorey | 
: cor { peo: | S96 | L66°  FI00 TOFOO 8000 F£SSIO | TIR | 2100 FLSEO | C'9F FOOD TFLIO SOF o een 

5: | | e 
: FCI0 =" GFIO =" oreo | 9620 ="F 

= SIL | ese | SLE | 866 | LO00 F 9900 | £000 FF600 | SES | 6000°+EEZ0" | ESS z000 FSET O0 | FEE SDEpD]y | 
$ see | ses: | CLE | 066 | stoovF+FZ8S00" | 6000 72600" | 609 | 6000 TOFD | 6°S9 £000 FR£IO : T09 mPa 
= izg—| 086 | LES | SEG | L000 FL I00 —| FOOD FLETO" LoL 0100 +6610" ESL 000 + 6II0 OL mv] 
= een | see | aLe | S66 | FIOO F 9800" | L000°+CSI0" TIF 2100 9680 | OTE LOO + ISIO | 9 CF EMBSEUrE X 
sso 196 | €LS | 968° LI00 T9600 OI00 FELIO «ges STOO FESQed | OTF  £000 TSEC | FFE _ ofgoI 
QI6 © 686 | 986 SSE | L 100 + 6610" S000 FCTIO | LE | CI00 F99€0 Tie F000 FEZE0 TLE IWUDTYD A 
gre | 298° | 296 | S86 | LS00 F5 900" s100 EF610 | att | 6100 F1SE0 | EF | LO00 7 160" | 905 | SEnoY uoreg 


LEQ” 098" 626 | 06680 Lz00°-7 LYTO FI00 +6FIO | IS | CIOO+E6S0° | FE C000 + L£EO | SOE SPuDS 2254 Al 
SOL es | 086 | 066 | TE00°+ 9610" 9100 T£SSIO Se . 100 TE9F0 96€ 9000 7 SSE0 SFE | SEM US 
ELL Soe : QS6 FEE GEO + COL" C100 TG66I0 TIE 9100 6SFO | GFE 9000 TOISO | 908 wophurysn 4 
2027 OFS OSE SSE 6 S00 + FEO" 0800 9610 CVE LI00 TO8C0 | 6G£E  S000 F68€0 | SEI FINQRIJ 
age” ESS" | £26 | S96%| 0 200°+0 10° IF0O-+A43FIO 0S ET 0200 + 19500 QE 9100 + Fresco | CLI ETES 
est Sai st a : g “a g “a 2 "N x 
SSA 7 Kg a, | == ; = se - £IOIRASYO 
CrFé+ ag + D'N=N (Fo+D'N=X (Ze + UN = 
| IA 


THE WORLDWIDE DISTRIBUTION OF THE F2 LAYER ELECTRON DENSITY ETC. 231 


tudes; however, this fact is not observed at Hobart and, on the other hand, 
it is not supported by the latitudinal variation of the partial correlation co- 
efficient 7,5... 

Concerning the latitudinal variation of « and 6, there appear some ana- 
logies, obviously less evident, with the variation of y, while the coefficients 8 
can be considered fairly symmetrical in the two hemispheres and little de- 
pendent upon the latitude. 

These features of «, 6, y, 6 can better be seen, if we assume the mean 
values &0, Bo, Yo, d a8 «normal » values and o,, Og; Os Os AS «normal » devia- 
tions: it clearly appears that the more significant dependence upon the lati- 
tude is that of the coefficient y, which is larger in the northern hemisphere 
and increases for increasing values of the geomagnetic latitude. 

In regard to the correlation coefficients, one sees that, again with the 
exception of Maui and Baton Rouge for the northern hemisphere, the values 
of 73. are clearly greater in the northern than in the southern hemisphere. 

The extended results of Table II seem to support the preliminary conclu- 
sions of a preceding paper (MARIANI (°)) of some different dependence or « sen- 
sibility » of the two hemispheres to the activity of hydrogen filaments (or to 
some other solar phenomenon closely connected). It may be interesting to 
compare our results with those of other authors: for example ALLEN (*°) who 
considered only six observatories (Fairbanks, Brisbane, San Juan, Maui, San 
Francisco, Baton Rouge) without any distinction of the northern from the 
southern ones, obtained for x a mean value «= 0.0200 which can be consi- 
dered in good agreement with our results (our mean value for the four obser- 
vatories of Brisbane, Maui, San Francisco, Baton Rouge gives «= 0.0211); 
with regard to the results of GALLET and RAWER (1), the values calculated 
by their formula giving the dependence upon the magnetic inclination are, 
respectively, for Freiburg, San Francisco, White Sands, Wakkanai, Tokyo and 
Yamagawa, the values «= 0.0186, 0.0185, 0.0183, 0.0178, 0.0161, 0.0156 
which do not agree well with our values. 

With regard to the values of N), the statistical errors of which are of the order 
of a few per cent, it appears that they are not rigorously symmetrical around 
the geomagnetic equator: Fig. 8 shows the values of N, for the case of mul- 
tiple correlation arranged, in the upper part according to the geomagnetie 
latitude (black dots) or the magnetic inclination (crosses), and in the lower 
part according to the sum of geomagnetic and geographical latitudes (black 
dots); in this latter case, a greater symmetry of N, is apparent; this fact is 
physically interesting because it shows that, in effect, although the abscissae 


(*) F. MARIANI: Journ. Atmosph. Terr. Phys., 10, 240 (1957). 
(10) C. W. ALLEN: Terr. Magn., 53, 433 (1948). 
(1) R. Garter and K. Rawer: Ann. Geophys., 6, 104 (1950). 
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are somewhat conventional, the electron density in the #2 layer varies either 
with the geomagnetic or with the geographical latitudes (i.e. with the mean 
zenithal angle); moreover, some of the 
coefficients «, y, 6 appear in a better 
order, if we rearrange the observato- 
ries according to the sum of the two 
latitudes. 

A singular case is represented by 
the values of N, at Christchurch which 
appear smaller than the others: this 
fact occurs systematically in all the 
single and multiple correlations, so that 
one is led to assume some particular 
local situation; the smaller value of N, 
accounts in some way for the greater 
values of a, f, y, 0. 

Concerning the latitudinal variation 
Fig. 8. — The dotted lines are symme- of N,, we can also consider the result 
trical lines approximately fitting the of BHAR (2) who has depicted the true 
calculated values. The northern lati- : 

A VA RT relation between noon values of foF2 
and the magnetic dip for the last sun- 
spot minimum, by eliminating the effect 

of solar zenithal distance variation: he plotted the noon values of foF2 against 
the magnetic dip for many observatories, keeping the solar zenithal distance y 
as a fixed parameter (instead of the equinox or any other time of the year); 
there remains an asymmetry in the two hemispheres which seems to increase 
for increasing values of y,,,, (particularly for y,,,.= 20° and = 30°): more 
precisely, there appears a maximum of foF2 at about 30°- 40° of magnetic 
dip, larger in the northern than in the southern hemisphere. 

Our results seem instead to suggest a more regular variation of N, exhi- 
biting two fairly symmetrical maxima. 

In order to infer some further information on the north-south asymmetries, 
the simple and multiple correlations of the monthly values of N—N,, with the 
monthly values of À and A, have been calculated. The seasonal and non- 
seasonal variations of N have been subtracted, because it seems that these 
variations, although a priori they may be correlated with the solar activity, 
are strongly influenced by factors other than the solar activity (i.e. the va- 
riation of the zenithal angle of the Sun, the height variation of the generalized 
recombination coefficient, etc.). In these new correlations, three linear re- 
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(12) J. N. BHAR: Journ. Atmosph. Terr. Phys., 10, 168 (1957). 
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gression laws have been assumed, expressed in the form: 


(4) MEN (NN al + «'R); 
(5) N Nig = (NN) (1 -+ 0’A,) ; 
(6) N — Nj, = (N — Nuh (1-+f'R+y'4,). 


The deductions one can draw from the results of the calculation, collected 
in Table III, substantially agree with those drawn from Table II; the simple 
and multiple correlation coefficients now seem fairly symmetrical in the two 
hemispheres; as regards the partial correlation coefficients, there appears some 
difference between the two hemispheres: Masi as a whole, is slightly greater in 
the southern hemisphere; 13, on the other hand, is symmetrical for lower 
and mean latitudes, while, for higher northern latitudes, it reaches values 
rather greater than at symmetrical southern latitudes. 

With regard to the coefficients «’, B’, y’, 0’ there appears again some clear 
tendency for «’, 6’ and particularly for y’ to increase in value for increasing - 
northern latitude; in this case, the values of y’ appear statistically significant, 
also in the southern hemisphere; moreover, one sees again the small depen- 
dence of £' on latitude and its similarity in the two hemispheres. 

With regard to the results relating to Christchurch, the values of (N — Ns)o 
are again smaller than the other values. 

The general diminution of the correlation coefficients is obviously due to 
the effect of the variables that we have excluded from our calculation; instead, 
a rather surprising feature is the remarkable increase of the values of the 
partial correlation coefficients 7,;,, with respect to the values of Table I, for 
the southern hemisphere. 

Finally, we have considered the correlations of the series of N — N,, —N 
amiga À, Al A ry) according to the linear regression laws: 


(7) NN, N = (NN, WN). + a(R — E); 
(8) . N—WN,,—WN = (N—N,,—WN).+ d(A,—A,) ; 
(9) N—N,,—WN = (N—N,,— WN), + 0(R — R)+ e(A,—A,). 


In these correlations, obviously, the values of (N —N,,—N), are very 
small and oscillate around zero. In Table IV we show the ratios a", 6”, y", 6” 
of the coefficients a, b, c, d to the corresponding values of (N — N»), obtained 
in the regression of N — N,, on R and A,: in this manner, the values of 
a, b, c, d appear more or less conveniently « weighted »; moreover, we have 
indicated in the Table the correlation coefficients. 
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Although the statistical errors on the ratios «", B", y”, 6” are very great, 
there remains however some indication of an asymmetrie « sensibility » of the 
two hemispheres to the hydrogen filaments. 


42. Second case: at midnight. — The study of the correlations can be carried 
out in a similar manner, using the corresponding midnight series of N. 

There again appears the discontinuity of the slope about the beginning 
of 1951: however, the dependence of N upon the areas A, is noticeable also 
in some observatories where it is fairly negligible at noon, so that, as a whole, 
the systematic asymmetries found at noon are practically absent; however, 
we must consider that, in effect, the midnight electron density depends upon 
the diurnal ionization intensity but, as in the phenomena of hysteresis, it is 
highly influenced by the foregoing behaviour of the diurnal F2 layer so that, 
from the physical point of view, the consideration of the correlations of nightly 
electron density with solar activity cannot be very significant. If, instead, 
we fix our attention on the behaviour of the maximum electron density, as 
a whole, we may assume that the asymmetric properties of the #2 layer in 
the two hemispheres are substantially and clearly present during daylight. 


4°3. Discussion. — From the foregoing analysis, it may be concluded that 
the dependence on the solar activity of the maximum electron density is fairly 
well expressed by the double correlation of N with À (or A,) and A4,; it is 
very interesting to notice that the north-south asymmetries of the correlations 
with solar activity, especially the one with the hydrogen filaments during 
daylight, show some interesting parallelism with the north-south asymmetries 
in the behaviour of N,,, which are present only during daylight. 

In order to explain the above mentioned asymmetries, the presence of some 
solar radiation could be assumed, principally or exclusively reaching the nor- 
thern hemisphere: if one excludes the direct effect of an electromagnetic ra- 
diation, which is more or less symmetrical around the Sun’s and the Earth’s 
equators, One could imagine some continuous asymmetrical corpuscular ra- 
diation; if, on the other hand, one excludes that this corpuscular radiation 
may reach the northern and the southern hemispheres with different intensity, 
as an effect produced by the geomagnetic field, one could imagine some se- 
lective effect determined, for example, by some geoelectrie field: in this case, 
considering for example a dipole field, one might be forced to assume a primary 
corpuscular radiation composed of only positive (or only negative) particles 
or, in every case, having a net (positive or negative) charge different from zero. 

Now, as it was pointed out by SIMPSON (*), if a variable geoelectric field 


(19) J. A. Simpson: Ann. Geophys., 11, 305 (1955). 
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is to be assumed, in order to explain the cosmic radiation variations, the neces- 
sary geoelectric potential is of the order of 2-107 to 4-108 V. A Coulomb field 
is excluded for, as FERRARO has shown, its potential cannot be higher than 
10? V, otherwise, with higher fields, the ionospheric layers would tend to blow 
away; there remains the other alternative of a radial electric field produced 
by ring currents or cylindrical current sheets, as considered by NAGASHIMA (14) 
and others. However, one has to remember that these calculations substan- 
tially refer to the disturbed conditions associated with magnetic storms and, 
on the other hand, that the calculated fields appear at least of an order of 
magnitude too small to account for the observed variations. Further arguments 
indicate the hypothesis of a geoelectric field as being highly improbable: indeed, 
the minimum energy necessary for a primary particle in order to reach the 
fairly low geomagnetic latitude of 40°--50° is of the order of some GeV; such 
a degree of energy is sufficient for the penetration of the atmosphere, so that 
the cosmic ray intensity at sea level, or, better, at high latitude, would exhibit 
some asymmetry in the two hemispheres, contrarily to the known experi- 
mental data. Moreover, many authors (FIROR (!), etc.) have shown that the 
geographical distribution of the intensity of cosmic rays, as it has been pos- 
sible to confirm on the occasion of solar flares, sufficiently agrees with the 
expected distribution of an extraradiation coming from a solar source and 
substantially subjected to the influence of only the geomagnetic field. 

An other possibility is the presence of some low energy radiation which is 
continuously accelerated near the Earth to the ionosphere by some local 
electric field and whose effects interest only the upper atmosphere where it 
is absorbed. In this case, some relation of this radiation with the aurora could 
be present. 

If this would not be the case, one could perhaps attribute, as in the case 
of the asymmetric features of N,,, the origin of the observed latitudinal ano- 
malies to some terrestrial effect, more indirectly influenced by the solar 
activity. 


5. — Further discussion and conclusions. 


As it has already been said, there appears some parallelism between the 
latitudinal asymmetries of N (or N — N, or N— N,,— N) and N,,. A de- 
pendence of N and N,, upon solar activity is indisputable; however, some 
remarkable features of the latitudinal dependence may seem influenced in a 


(4) K. NagasHIMA: Journ. Geomagn. Geoelectr., 3, 100 (1951). 
(15) J. W. Frror: Phys. Rev., 94, 1017 (1954). 
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rather indirect manner. A tentative suggestion is that of considering the effects 
induced in the ionospheric electron densities by the general circulation of the 
upper atmosphere. 

We, at first, consider the effects of the sphere tides: enough obser- 
vational evidence on tidal phenomena is furnished by current literature, espe- 
cially with regard to the diurnal anomalous effects on foF2. Martyn’s theory 
(MARTYN (1f)) of atmospheric tides in the ionosphere can give a solid base 
for accounting for such anomalies in à determined place; however, if one ac- 
cepts his point of view, in order to explain the latitudinal variations, and 
particularly the north-south asymmetries of N,,, one has to imagine some 
asymmetric latitudinal characteristics of the tidal movements. As MARTYN (1) 
has shown, in order to explain the behaviour of foF2 at Cape York, Brisbane, 
Canberra and Hobart, one can assume the presence of two spherical harmo- 
nics 43, y; (semidiurnal components) in the general expression of the velocity 
potential. At first sight, one might think that the features of N,., particularly 
those relating to its phase-angle, varying so differently with the latitude at 
noon and at midnight, could be explained by assuming also the presence of 
the diurnal terms y, y, in the velocity potential considered above. 

We have thus attempted to frame in some detail the experimental behaviour 
of N,, in such a scheme, assuming the presence of the terms wi, Py, Wo, Ys 
in the velocity potential; however, we have not been able to obtain some 
complete agreement of experimental and theoretical results. In effect, it must 
be considered that the terms »? present in the general expression of the velo- 
city potential are seasonal or annual terms, while, as a result of the discussion, 
of Section 3, this does not seem the case of N,,. On the other hand, the cor- 
rections introduced by the presence of the geomagnetic field render the tidal 
circulation substantially symmetrical around the geomagnetic instead of around 
the geographical equator, because the effective asymmetries they introduce 
in the two hemispheres are negligible. 

There remains the possibility that a mechanism of the Martyn tipe is 
effective, but that the circulation in the upper atmosphere is not well expressed 
in terms of spherical harmonies potential or that an asymmetric, non- Une 
conductivity is present. 

In other words, it may be suggested that the observed north-south asym- 
metries of N,,, may be interpreted as the effect induced by the geomagnetic 
field on the worldwide circulation of the upper atmosphere, which may be more 
‘or less different in the two hemispheres. In effect, the «regional» behaviour 
we have observed in N,, may perhaps be an indication of some regional cha- 


(18) D. F. Martyn: Proc. Roy. Soc., A 189, 241 (1947). 
(7) D. F. MARTYN: Proc. Roy. Soc., A 194, 445 (1948). 
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racteristic of this general circulation. Moreover, other experimental results 
seem to support the hypothesis of some different characteristic of the northern 
with respect to the southern hemisphere: on discussing the anomalous varia- 
tions of the maximum electron density in the F2 layer, SATO (#8) also has 
indicated some asymmetries between the two hemispheres in the ratio > of noon 
median values of the electron density in the minimum cos y month to that in 
the maximum cos y month, at various latitudes. This ratio is greater than unity 
for northern latitudes higher than 10°, at sunspot minimum, and becomes 
larger for higher northern latitudes, at sunspot maximum. The ratio r is larger 
for higher geographical latitudes, more for northern than for southern lati- 
tudes; moreover, the seasonal anomaly is amplified on a worlwide scale during 
sunspot maximum far more than during sunspot minimum. SATO has finally 
shown that the seasonal anomaly of foF2 cannot be accounted for merely by 
the vertical drift. 

On the other hand, SCHWERDTFEGER and PROHASKA ("*) have pointed out 
the occurrence of remarkable latitudinal asymmetries on the mean of the 
amplitude of the semiannual pressure oscillation and of its ratio to the ampli- 
tude of the annual oscillation. 

Obviously, the above outlined point of view requires further observational 
data. The most important questions which must be exhaustively studied are: 


i) the behaviour of N and N,, and the correlations with solar activity, 
for observatories located in the equatorial belt and at high northern, and parti- 


. cularly at high southern, latitudes; 


ii) the worldwide distribution of the ionospheric winds, their diurnal 
and seasonal or non-seasonal variations; with regard to this question, some 
other research on the regular (not necessarily seasonal) tidal movements can 
be very useful. Moreover, one cannot disregard the possibility of a more or 
less regular slipping back of the upper atmosphere, with respect to the motion 
of the lower atmosphere and of the Earth: this effect should equally interest 
both hemispheres. 


It will thus be very advantageous that the intensified researches on iono- 
spheric physics (and on allied branches) during the International Geophysical 
Year may be continued regularly in the future. 

An extensive statistical study of all the ionospheric data for the years 
1938-1958 and of their correlation with other significant solar parameters is 
in preparation. The numerical calculation is now in progress by the Electronic 
Computer of the « Centro di Studi sulle Calcolatrici Elettroniche » of the Uni- 
versity of Pisa. 


(18) T. Saro: Journ. Geomagn. Geoelectr., 6, 99 (1954). 
(1°) W. SCHWERDTFEGER and F. ProwasKa: Journ. Met., 13, 217 (1956). 
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RIASSUNTO 


L’andamento globale della variazione di periodo 12 mesi della densità elettronica 
massima nello strato F2 risulta sensibilmente differente a mezzogiorno, e a mezzanotte. 
In ogni caso, essa è controllata dal campo magnetico terrestre e dall’attività solare 
ma, mentre a mezzogiorno mostra notevoli asimmetrie tra i due emisferi, al contrario 
a mezzanotte, le sue caratteristiche sono simmetriche. Non risulta possibile attribuire 
tale diversità di comportamento della parte illuminata e della parte oscura della Terra 
all’effetto di qualche sorgente esterna di radiazione. Si è indotti a pensare in termini 
di qualche asimmetria di origine terrestre: per es. di asimmetrie nella circolazione 
generale dell’alta atmosfera. Quanto alla correlazione con l’attività solare, di nuovo 
appaiono sensibili asimmetrie tra i due emisferi N e S; in particolare, l’emisfero N 
appare maggiormente soggetto all’influenza dei filamenti cromosferici solari, in misura 
tanto maggiore quanto maggiore è la latitudine geomagnetica. L’interpretazione di 
questa asimmetria in termini di una diretta influenza del Sole sullo strato F2 costi- 
tuisce un grosso problema; si è indotti a pensare a qualche effetto più o meno scono- 
sciuto di una radiazione corpuscolare (o elettromagnetica) incidente sull’alta atmosfera, 
che ne influenzi per esempio i movimenti o la conducibilità. Allo stato attuale delle cono- 
scenze, in attesa di ulteriore e qualificato materiale sperimentale, non si può escludere 
la possibilità che ambedue i tipi di asimmetria, quella relativa alla variazione di pe- 
riodo 12 mesi e quella relativa alla correlazione con l’attività solare, siano determinate 
da una medesima causa. 
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Summary. — Interactions of K*-mesons with emulsion nuclei have been 
observed in a scan-of 185 m'of track. The analysis yielded information 
on scattering reactions K*-p and K*-n and charge-exchange reactions. 
An energy dependence of the three cross-sections up to 350 MeV was 
obtained, from the collected data of different laboratories working with 
emulsion technique. The phase-shift analysis in the 300 MeV region 
shows a contribution of 7=0 state of the same order of magnitude as 
that of T7=1 state. An interference seems to be present between the 
P-waves of the two states, depressing the K*-n scattering cross-section 
and raising the charge-exchange cross-section. 


1. — Introduction. 


The results of many investigations on the interactions of K*-mesons in 
nuclear emulsions in the region from very low energies up to about 200 MeV 


n 
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have been published recently (1-1). The results have been found to be well 
described in terms of the optical model for complex nuclei. From the inter- | 
actions with the Hydrogen and other nuclei in the emulsion, certain features — 
of the elementary interactions with free protons and «free» neutrons were 
inferred notably the slow rise in the K* proton interaction (7 = 1 only), 
and the much faster rise in the K* neutron cross-section (T — 1 and T= 0) 
from a value at low energies which suggested that there was probably only 
a very small 7= 0 part, 

The present experiment was designed to extend, by an essentially similar. 
technique, the investigation of the K* meson interactions to the energy region 
(240--300) MeV. Attention was again concentrated on the magnitude of the 
interaction cross-section as a function of the energy, the ratio of the number 
of Charge-Exchange to non-Charge-Exchange events, and? the collisions with 
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E. ILorr, A. PEVSNER and R. Rirson: Phys. Rev., 101, 1617 (1956). 

(2) N. N. BIiswas, L. CECCARELLI-FABBRICHESI, M. CECCARELLI, K. GOTTSTEIN, 
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(5) B. BHOWMIK, D. Evans, 8. Nirsson, D. J. PRowsE, F. ANDERSON, D. KEEFE, 
A. KERNAN and J. Losty: Nuovo Cimento, 6, 440 (1957). 

(6) M. Wipcorr, A. PEVSNER, D. Fo uRNET-DAVIS, D. M. Rirson, R. SCHLUTER 
and V. P. Henry: Phys. Rev., 107, 1430 (1957). 

(?) T. F. Hoane, M. F. Kapton and R. CestER: Phys. Rev., 107, 1698 (1957). 

(8) D. FouRNET-DAVIS: Phys. Rev., 106, 816 (1957). 

(®) C. MARCHI, G. QUARENI, A. VIGNUDELLI, G. D’Ascora and S. Mora: Nuovo 
Cimento, 6, 1790 (1957). 

(1°) B. SecHI-Zorn and G. T. Zorn: Phys. Rev., 108, 1098 (1957). 

(11) J..M. LANNUTTI, S. GOLDHABER, G. GOLDHABER, W. W. CHUPP, S. GIAMBUZZI, 
C. MarcHI, G. QUARENI and A. WATAGHIN: Reports of the Padua-Venice Conference 
(September 1957), III, 1. 

(2) M. GrILLI, L. GUERRIERO, M. MERLIN and G. A. SALANDIN: Reports of the 
Padua-Venice Conference (September 1957), III, 16. 

(1*) N. N. Biswas, M. CeccARELLI and N. ScHMITZ: Reports of the Padua-Venice 
Conference (September 1957), III, 25. 

(14) L. T. KeRTH, T. F. KyciA and L. Van Rossum: Reports of the Padua- Venice 
Conference (September 1957), III, 28. È 

(5) G. Ico, D. G. RAVENHALL, J. J. TIEMANN, W. W. CHuPP, G. GOLDHABER, 
S. GOLDHABER, J. M. LANNUTTI, and R. M. THALER: Phys. Rev., 109, 2133 (1958). 

(16) B. BHOWMIK, D. Evans, 8. Nirsson, D. J. PRowsE, F. ANDERSON, D. KEEFE, 
A. KERNAN, N. N. Biswas, M. CECCARELLI, P. WALOSCHEK, J. E. HoopER, M. GRILLI 
and L. GUERRIERO: Nuovo Cimento, 5, 994 (9157). 

(1?) M. GRILLI, L. GUERRIERO, M. MERLIN, and G. A. SALANDIN: Nuovo Cimento” 
10, 205 (1958). 


1154 


THE INTERACTIONS OF POSITIVE K-MESONS WITH NUCLEI ETC. 243 


Hydrogen nuclei in the emulsion. In addition, a new feature present at this 
energy was the possibility of x-meson production, the threshold for this process 
in the collision of a K-meson with a stationary nucleon being 220 MeV. Some 


| preliminary reports of this work have been given at the Geneva Conference 


(1958) and elsewhere (18-20), 


2. — Experimental details. 


21. Exposure conditions. — The plates used were from a stack (Kj) of 
240 Ilford G-5 600 um strips 20 em X17.5 cm, exposed to the Berkeley sepa- 
rated K*-meson beam (nominal momentum — 625 MeV/c), processed at Bristol, 
and shared among the laboratories of Bristol, Dublin U.C., and Padua. The 
number of K-particles entering each emulsion strip at the shorter edge was 
about 100. The x/K ratio varied from place to place in the beam and gene- 
rally was between 1:1 and 3:1. The development was reasonably high—a 
minimum blob-density of about 22.1 blobs/100 um—the development gradient 
was only a few percent at most, and the distortion was negligible for the range 
of energy determination considered in this experiment. In addition, the values 
of minimum blob-density in different plates were remarkably similar. 


22. Selection and analysis of events. — The tracks were chosen for following 
by scanning in a line at right angles to the beam direction. This line was 
chosen at 10 mm from the beam entry point for one part of the experiment 
and 15 mm for the other, in order to allow the identification of the primaries 
of stars found close to the pick-up line by multiple scattering measurements. 
Each track was blob-counted under X100 objective over a distance of 2 mm, 
which was sufficient for an efficient selection of the heavy meson tracks. The 
blob density at this momentum was very close to minimum (22.1 blobs per 
100 um) for K-mesons. As an additional check on the choice of track and 
the correctness of tracing through, a large sample of tracks was again blob 
counted after it had been followed through for about 8 cm where the expected 
ionization for K-particle tracks was about 31 blobs/100 um. 

In all cases, when an interaction was found, identification of the primary 
was demanded. If one of the tracks emerging from the star was definitely 
identified as a K-meson, the primary, must also have been a K*-meson. In 


(18) M. GrILLI, L. GUERRIERO, M. MERLIN, Z. O’FRIeL and F. A. SALANDIN: Nuovo 
Cimento, 10, 163 (1958). 

(19) D. KrEFE, A. Kernan and A. Montwitt: Nuovo Cimento, 10, 538 (1958). 

(2) D. Evans, F. Hassan, K. K. NaGPAUL, M. D. SHAFI, E. HeLmy, J. H. MULVEY 
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all other cases, scattering measurements were carried out on the primary 
track to determine its mass. Only in 8 cases out of 365 was the identity of 
the primary found to be other than a K-meson (5 x-mesons, 2 protons, 1 
cosmic-ray star with outgoing prong in the K-beam direction). To correct for 
tracks followed, which were due to other particles, the total track length was 
reduced by 2%. (The mean free path for 7-meson and proton interactions 
is somewhat smaller than for the K*-meson). 

In the case of each star found, every prong was, if possible, followed to 
rest and its end examined for evidence of decay. The energy and direction 
of emission of the prong was recorded. In cases where a particle left the stack, 
its mass was determined by measurement of increase of ionization as a function 
of range, or, if flat enough, by measurements of (%, b*). 


2°3. Estimation of the primary energy. — The mean energy of the incoming 
beam and its inherent spread in energy were estimated in two ways. Firstly, 
a sample of 45 K-meson tracks was chosen at random and they were followed 
for 20 cm—that is, nearly to the edge of the stack—and ionization measure- 
ments made. In this way, 

the mean residual range at 

that point was found to be 


12 ange adiciribution 9.2 em and the — 
10 for not interacting K's deviation of the range distri- 

bution 2.8 cm (see Fig. 1). 
i R=292228 mm Taking into account the er- 
i rors due to straggling and 
ia theionization measurements, 
2 the energy at entry of the 
9% 540 280 320 360 400 beam turns out to be 310 MeV 


and the spread — 20 MeV. 
The second method involved 
the estimation of the primary 
K-meson energy in each of the K-Hydrogen events by simply summing 
the energies of the outgoing proton and K-meson. The K-meson energy 
at the point of interaction can be corrected up quite reliably to give the 
energy at the entry point and each individual estimate is good to a few MeV. 
The distribution in values thus found for our 13 events is shown in Fig. 2. 
The mean value at entry derived therefrom was (303 + 4.5) MeV, and the 
spread, — 16 MeV. Accordingly the mean energy at entry was assumed to 
be 306 MeV. 

Since it was desirable to divide the data into two energy groups, £ 270 MeV, 
individual energy estimates were required for the primaries of every star and 
here somewhat different methods were used by the two groups. The Dublin 
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group took simply the expected value at the star, based on the distance of 
the event from the edge of the stack and the assumed average beam energy 
at entry. Thus each primary estimate is uncertain by about + 18 MeV. 


Entrance energy 
of K-H events 


È =(303+4.5)-MeV 


4 

O (E) = 16.3 MeV 
3 
2 
Pejo 
DEGREE 
200 220 240 260 280 300 320 340 366 


Fig. 2 


The Padua group determined the energy of the primary K-meson at each 
interaction directly by a grain-count along 2mm of track combined, when 
necessary with multiple scattering measurements over 2cm. The accuracy 
of this energy determination was + 6%. 

The energy distribution so determined for interacting K’s turns out to be 
in agreement with that deduced for non-interacting K’s, as was expected as 
the energy dependence of the interaction cross-section is slight. The experi- 
mental spread of + 6% coupled with the real energy spread of the beam 
results in some energy estimates less than 240 MeV and greater than 300 MeV. 
Because of such uncertainties it is clear that when a division of the data is 
made at 270 MeV, there will be an inevitable « spill-over » of events close to 
270 MeV from the lower energy category to the higher one and vice versa. The 
number of events involved is small and the effect on the results is negligible, 
since the mean free path varies only slightly, and the energy distribution is 
close to symmetric. 


3. — Experimental results. 


3 1. General considerations. — In all, 185 m of K-meson track have been 
scanned, and 13 K-Hydrogen events, 203 nuclear inelastic scatters, 91 charge- 
exchange events, 7 «stops», 37 decays in flight and 12 further inelastic un- 
analysable events were observed. The full details of these events are tabulated 
in Appendix I and the main features are summarized in Table I-a for the two 
energy intervals (240 270) MeV and (270--300) MeV. A «stop » is defined as 
an event where the K-particle stops without any prongs or decay product 
being apparent. It is reckoned that two of the observed «stops » were due 
to a decay-in-flight where the secondary particle escaped detection. The re- 
sults on K-H events are given in Table I-b. 


16 - Jl Nuovo Cimento. 
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From Table I it can be seen that the measured mean free paths for the 


various processes in emulsion are: 


M.F.P. 
M.F.P. 


for K-H collisions . 


for inelastie collisions 


(excluding K-H events) 
(scatterings + charge-exchange) 


14.2 GIS 


0.61 + 0.05 
0.58 + 0.05 


| 
| 
| 0.84 + 0.07 
“I 
| 
| 


(240 + 300) MeV 


(240 +270) MeV 
(270 +300) MeV 


40270) MeV 


) 
“ 
27 
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(*) Two stops were due to a decay in flight where the secondary particle escaped detection. 


TaBLE I (b). — K- Svago collisione. 


(: 
M.F.P. for inelastic scatterings . 
i 0.92 + 0.08 (270-300) MeV 
2.25 + 0.33 (240-270) MeV 
M.F.P. for charge-exchange collisions 
E 1.56 + 0.20 (270-300) MeV 
TABLE I (a). 
Energy interval | | 
i | LATI MINARE] Spurious 
pitts Cee teal cna mL à Events | 
(240-270) MeV (270-300) MeV | 
5 ad Baan aie aks lads (Adie fas nee | 
= = 2 | Ss 3 EX 2 a 
.. |SB/ Eg! ©| Ses) & |Sg|Egl © seg SSA 
(S| gf 3 | > | SD go | 388) .d 2 | 3 | 83 | SE | 
Labt.| $3|48| &|K-w| @|88| SE |SS 48) è KH Sa | SE ag | 
SalSh à s |Ae| 8° ASS) à s | Aa | à SE 
© D | Ay O 5 | By 
| Di = ina = ee = 

Pa. | 58 | 21 | 12 | 5 | 2 | 10 | 53.5 | 50 | 30 | — | 4 | 6 | 11 | 55.5 3 

Du 45218002 2 2 6 | 40.7-| 43 | 22 411182 le 10 [13670 5 
Total|110|39 | 3 7 4 | 16 | 94.2 | 93 | 52 4/6] 8 21 | 90.5 8 | 
(2) (3) | 


N tim het Number x 

of events E, Xcm of events Ei XCM 
Pd. 99 250 38° Pd. 678 230 76° 
Pd. 192 260 94° Pd. 1428 240 86° 
Pd. 99 250 38° Pd. 678 230 76° 
Pd. 192 260 94° Pd. 1428 240 86° 
Pd. 386 270 106° Du. 220/14 245 80° 
Pd. 620 232 135° Du. 181/7 274 108° 
Pd. 637 265 76° Du. 202/27 283 136° 
Pd. 802 325 44° Du. 217/21 233 98° 
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The Ratio of Charge Exchange/Non-Charge Exchange events in the lower 
energy is: 0.37 + 0.07 and in the higher energy interval 0.59 + 0.10. 


Only 


one case of z-meson production in a K*-nucleus collision was ob- 


served; this star had one prong—the emergent x*-meson—and the experi- 
mental details are described fully in (18). 

Fig. 3 shows the world data or the K-H angular distribution, dap/dQ, 
including the 13 events found in the , esent work. Fig. 4 is a scatter dia- 


gram of 
G. > for 


the inelasticity, AE/E, against the laboratory angle of scattering, 
the 203 inelastic events. 


K-H events 
200 < E, <350 


0 - -2A - cm. 
Fig. 3. 


4. — Discussion of results. 


~ 


41. The mean free path for inelastic scattering. — In Fig. 5 are shown the 


values at different energies of the mean free path of K* mesons in emulsion 
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obtained in the present work and at vàrious laboratories (K-Hydrogen events 
are excluded) (12°). The value of the mean free path is seen to decrease 
with energy, and at low energy (< 130 MeV) this decrease is well explained 
in terms of the Pauli exclusion principle which prohibits low momentum trans- 
fers to ‘bound nucleons. The continued decrease at higher energies is not 
explicable in the same terms but.is due as discussed below to an increase in 
the K*-bound-nucleon cross-section. 


42. The frequency of charge exchange interactions. — Fig. 6 shows the 
available world data on the ratio of charge-exchange (CE) to non-charge- 
exchange (NCE) interactions. The points at 245 MeV and 280 MeV have 


SI 
.. _ charge exchange 
80 Ratio = scattering 
| * UCLA 

70 © BNL 

o Pd 
60+ e Du+Br 

4 G6 
50r x Pd* Du 
40} 
30! 
20! 
10 
% 50 100 160 200 20 300 380 


been obtained in the present experiment. It can be seen that while all the 
results below 200 MeV could have been consistent with the value 0.2 predicted 
by the Charge Independence Hypothesis for the case of a vanishing T = 0 
interaction, the higher energy points lie substantially above this value and the 
rapid rise seems established. Since an emerging K-particle may have had more 
than one collision within the nucleus its chance of suffering charge exchange 
is greater than in just a single collision with an « average » bound nucleon. 
Fig. 7 shows the experimental curve of Fig. 6 after correction for double scat- 
tering—and represents the CE/NCE ratio for a single collision with an average 
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bound nucleon. At the energies considered in the present experiment double 
scattering occurs in about 40% of cases (21). 


._ C. EX 
=> for 
Ratio sai corected fo 


double scatterings 


Fig. 7. 


43. The K*-nucleon cross section. — Adding our data to those reported at 
the Geneva Conference the K-proton cross-section results to be (21.344) mb in 
the energy interval (200--350) MeV Fig. 8 shows the results on the total cross- 


150 
Fig. 8. 


section, o,, at different energies. It is likely that there is a genuine rise in o, 
in the high energy region. 

One can now deduce the K-neutron cross-section, o,, by using the mean 
free paths given above in Section 3. Firstly the M.F.P. must be decreased 


\ . (1) K. A. BRUECKNER, R. SERBER and K. M. Watson: Phys. Rev., 84, 258 (1951). 
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by a factor (1 —V,/E,) ~ 0.96, to correct for the effect of electrostatic re- 
pulsion, where V, = average Coulomb potential at the nuclear surface. 

In order to deduce the cross-section on the average free nucleon, allowance 
must be made for the effects of the Pauli exclusion principle and of shading. 
The former correction at these energies is less than 10% if the angular cross- 
section in the C.M. system is assumed to be not too far from isotropic but the 
shading correction is much larger and depends on the assumptions made about 
the nuclear size and shape. Estimates of nuclear sizes and shapes vary very 
widely and clearly have a precise meaning only in relation to the experimental 
context. Whether one takes a one-parameter or two-parameter «shape », it 
appears that the magnitudes of the parameters will depend on the energy and 
nature of the bombarding particle. Since there is no evidence for a long 
range K-nucleon force and the energy involved is fairly high, we have chosen 
the simple homogeneous model with radius R—=7,A*, and have calculated 
results for » =1.36 fermi, 1.25 fermi and 1.15 fermi. The middle value 
1.25 fermi, is in fair agreement with the high energy neutron and proton 
scattering data where a homogeneous matter distribution seems to be a good 
choice (cf. GLASSGOLD, ref. (?)). In addition the calculation has been carried 

out for a nuclear density 


i distribution: 

mb Cross sections for 
267 —K-bound nucleon sn AIA; Lo 

---K-free  , se te er) | 1+ exp(R aie 

oul + exp [ o/d] 
0 with R,—1.044#fermi and 

O RON PR ee ne 6425/ d= 0.57 fermi as used by 
161 IGo et al. (5). The results 


12 ee EA A I i 6=137f based on this shape factor 
are very close to those cor- 


8r responding to r,—1.25 fermi. 
n Fig. 9 shows the K*-nu- 
5 MeV cleon cross-section deduced 


RI, 
0 50 100 150 200 250 300 350 a from data reported in Fig. 5 
as described above for diffe- 
rent assumptions and one 
can see how sensitive the 
magnitudes are to changes in 7). The curves show the variation in o(N) 
before and after correction for the effects of the Pauli exclusion principle 
(isotropic cross-section). For the following calculations we have assumed 
% = 1.25 fermi. 


Fig. 9. 


(22) Intern. Congress on Nuclear Sizes and Density Distribution, in Rev. Mod. Phys., 
30, April (1958). 
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Using now the ratio (Charge-Exchange/Non-Charge-Exchange) at different 


energies, 
scattering—which can occur 
from either protons or neutrons 
—and the other corresponding 
to charge-exchange—which can 
of course arise only in a col- 
lision with a neutron (see 
Fig. 10). So by this means 
we can obtain directly o,(CE), 
the partial cross-section on a 
neutron for charge-exchange 
scattering alone. 

In order to determine the 
other partial cross-section on 


neutrons, that is, for non-charge exchange scattering, 


o(N) has been resolved into two parts, one corresponding to inelastic 


ox bound IA RT sr 


20+ o,-free 
mb} 0x-bound 
16+ ay-free 


12- 630125 10° cm) 


(ch ae 


” » oe os 


8 
4 sat eee 
0 one re Pa L i Me E 
0 50 100 150 200 250 300 350 
Fig. 10. 
o, (NOE), one can use 


the fact that the atomic nuclei in emulsion are composed approximately of 


54% neutrons and 46% protons. 


Thus 


o(N) = 0.460, + 0.54 6,, 


where 


On = CAC pe 


E) = Total cross-section on a neutron. 


Using the cross-section on an average nucleon corrected for the Pauli prin- 


ciple, o(N), and the K-p cross-section from the Hydrogen scatters, o,, 


ox” bound neutron (scatt) —— 
a, free neutron (scatt) ---- 


16 € (A =125-10 "cm 


(9, const) 
(Fy, pincreasing) 


Fig. 11. 


(b) 


that 0, 
15.7 mb. The behaviour of o,(NCE) so deduced is shown in Fig. 11 and 
it clearly passes through a maximum around 220 MeV independently of the 


n 


150 200 250 300 350MeV” 


the 
total cross-section on a «free » 
neutron, o,, may be deduced. 
Thus o,(NCE) may also be 
derived. 

However, the data of 
Fig. 8 cannot be taken as 
positive proof of a rise in 
o, with energy, so the cal- 
culation of o,(NCE) has been 
carried out under the two 
assumptions: (a) that o, in- 
creases with £, according to 
the data quoted in Fig. 8; 


is independent of energy below 300 MeV and has the value 
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assumption (a) or (b). At present it appears that assumption (a) is a bet- 
ter approximation to the truth than (b). 


44. The K-nucleon differential cross-section. — The distribution in center- 
of-mass scattering angle y.y has been shown in Fig. 3 for the 13 K-Hydrogen 
collisions reported in the present work and 13 other events found in other 
laboratories. The forward/backward ratio is 3 for the Du-Pd data and 10/16 for 
the world data. 

In the high-energy inelastic interactions of K*-mesons with nuclei the 
picture of collisions occurring with just single nucleons within nuclear matter 
seems to hold quite well so that in principle one could learn something of the 

differential cross-section on the «average » bound 
— Exp. distrib. nucleon. Apart, however, from the complication 
Zor -—- Calc. distib. that a nucleon is in general moving with a consi- 
derable momentum when struck, in practice a 
substantial fraction of the K*-mesons undergoes 
more than one collision within the nucleus. In 
Fig. 4 it can be seen that over 30% of the points 
lie outside the limits of (AE/E, 0) allowed for 
collisions with a nucleon of a Fermi gas (Pra 
= 241 MeV/c) and presumably they correspond 
to multiple collisions. The occurrence of a large 
amount of secondary scattering will smear consi- 
derably the laboratory distribution of single scatterings and this effect is 
difficult to estimate. The distribution in the values of 6 observed is shown 
in Fig. 12 together with that computed for single scattering from a Fermi 
gas, assuming (do/dQ),,, to be isotropic. The calculated distribution is slightly 
more peaked than the experimental one, but when the smearing effect due to 
multiple scattering is taken into account the calculated distribution will be 
less peaked, and that seems to indicate a (do/dQ),,, (K-free nucleon) a 
litte peaked in the forward direction. 


45. Decays in flight. — 37 examples of K*-mesons decaying in flight with a 
visible secondary track have been observed. From the depth distribution of both 
the decays in flight and the stops it is estimated that two events classified — 
as «stops » were decays in flight in which the secondary product escaped ob- 
servation. 

Thus the estimated total number of events occurring in 185m of track 
is 39. Since the mean velocity of the K mesons is 0.764 c, the lifetime turns 
out to be (2.1 + 0.3)-10-§ s seen in the laboratory frame. Dividing this by 
the time dilatation factor 1.55 gives the lifetime estimate for K* mesons at 
rest (1.37 + 0.22)-10-° s. 
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46. A tentative phase shift-analysis. — According to the Gell-Mann—Nishi- 
jima scheme, the only reactions expected from conservation of isotopic spin 
and strangeness considerations, are the following: 


K*+p->K*+p T=1 
er Sens eee T=1 and T=0 

K'+n+K°+p T=1 and T=0. 
If we assume that only S and P, and P; waves are responsible for the 


interactions and define 
Goo = EXP [4090] SID doo 3 


a the scattering amplitude in the 7 = 0, 2=0, state (the second indices 0, 1, 3, 
lignify S, Pi, Pi waves respectively), the total and differential cross-sections 


‘ san be written in the following way: 


(A) o = 4722(A 4-40), 
(B) a = À(4 + B cos 0 + C cos? ?), 


where 224= de Broglie wave length and 


(i)) A, = (jo)? + (@ig — M1)? 
(i) B_—2-Re[a}, (20,3 + au)] 
(iii) C3 = (248 + Au) — (Q13 — Au)? 


for proton scattering and 


(iv) Ale = (G10 + doo)? + À [ (ais — a) + (Gos— 1) |? 

(v) Bi = $ Re (ayy + doo)" [ (213 + dx) + (os + do1)] 

(vi) Ce = [Cas + du) + (2@0s + 1) |? =F [ (ais — @1) + (dos — 1) |? 

for neutron scattering, and neutron charge-exchange respectively. (Positive 
signs for scattering). 


Low energy results. From low energy data obtained in previous 
work (45161728) we can assume that: 


(2) O. R. PRICE, D. M. Srork and H. K. Trorro: Phys. Rev., Letters, 1, 212 (1958). 
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1) The most important contribution to the K-nucleon cross-section is 
from the T— 1 state in the S-wave. 6,,< 0 because the nuclear potential 
is repulsive. + 


2) The contribution from the T—0 state is smaller and is from a 
P-wave which gives a backward peaked K-neutron cross-section. Therefore 
do, must be positive (Eq. (v), P, terms neglected). 


High energy results. In analysing the high energy data, we can tentatively 
assume that the contribution of the P; waves is negligible. With this hypo- 
thesis, we can justify the maximum found for the K-neutron scattering cross- 
section at high energy, by assuming that the P-wave in the T— 1 state 
increases rapidly for energies > 200 MeV. 

If the phaseshift 6,, is negative, destructive interference will result bet- 
ween 4; and a,, which will lower the K-neutron cross-section for scattering. 
Constructive interference will take place for charge exchange (change of sign 
in (iv), (v) and (vi)). This agrees with the experimental data (Fig. 10). 

The experimentally found increase in the K-H cross-section supports the 
hypothesis of a rapid rise of the a,, contribution for E, greater than 200 MeV. 
Other supporting evidence for the above is obtained from the forward peak 
in the differential cross-section for K-nucleon scattering in the laboratory 
system (Fig. 12). This distribution was flat at low energies. 

Assuming that 6,,=—64,, in the energy region ~ 300 MeV, (as the neutron 
scattering cross-section seems to have a very small value (Fig. 11)), it is pos- 
sible, from the charge-exchange cross-section at this energy (Fig. 10) and cor- 
recting for the Pauli principle, to obtain a value for the scattering amplitude. 
Thus we obtained a P-wave contribution for proton scattering at 270 MeV 
of — 6.3 mb (taking 7) =1.25-10-!% em). On adding this to 13.5 mb, the as- 
sumed S-wave contribution, a value of 19.8 mb obtained for the expected 
K-H cross-section at this energy, in good agreement with the experimental 
value of (21.3 + 4) mb. 

Using the K-H cross section results, with an S-wave contribution of 
~ 13.5 mb, it is possible to calculate the P-wave phase-shift, in order to obtain 
a cross-section of 21 at high energy. We have thus calculated the expected 
K-H differential cross-section in the C.M. system. This gave a pronounced for- 
ward-peak, which does not agree with experimental data at present available (*). 

Similar conclusions can be obtained by considering the P, contribution 
to be negligible in comparison with P,. The main difference lies in the fact 


(*) Note added in proofs. — In a recent paper (74), indications are ise in favour 
of a forward peaking in the K~-p cross-section. 

(4) T. F. Kycia, L. T. KeERTH and R. G. BAENDER: Bull. Am. Phys. Soc., 4, 

5 (1959). 
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that term © in eqn. (iii) no longer vanishes, and an even more anisotropic 
K-H differential cross-section is obtained. 

In addition to the above difficulty, we have had to assume that the P-wave 
in the T=1 state, with a negative phase-shift (0,3) corresponding to a repul- 
sive potential, increases more rapidly than the P-wave of the T— 0 state, 
which has a positive phase shift. This point is difficult to interpret theore- 
tically. 


5. — Elastie scattering. 


In a sub-sample of 75.7 metres of track followed, a total of 167 elastic 
scatters was observed with plane projected angles greater than 2°. These 
were classified in intervals of the space scattering angle, and the differential 
cross-section derived for each interval. 


à | i | » 

| Redi ce 7°= | ge |10°+ |19°+ |14° : 

= Leo 7e) <8" | da 12°| -14°|=-19° 

© ee 
No. of events LB RIA 22 ISIS AT. 13 [17 | 20 BAR Hie 

| 10-* do/d® (mb/sr) | 59.5 | 14.3 | 10.3 | 6.8 | 450232 LT Lo) 0:6, 0.4 


Using the optical model, and assuming various values of the nuclear po- 
tential, theoretical curves of the elastic differential cross-section have been 
calculated by Evans et al. (2°). The points obtained in this experiment lie 
on an average about 25% higher than those found by these authors in a scan 
of 49 m of track and therefore reinforce their argument in favour of a sub- 
stantial repulsive potential at these energies. 


6. — Conclusions. 


1) At higher energies the contribution of the 7 =0 state becomes ap- 
proximately of the same importance as that of the 7’= 1 state. This fact can 
be deduced from the rise in the charge-exchange cross-section which reaches 
approximately the same value as the K-H cross-section (Figs. 8 and 10). 


2) The K-H cross-section seems to increase as the K-meson energy 
increases (from ~ 14 mb at (20 +200) MeV to ~ 21 mb at (200--350) MeV). The 
errors in the determination are rather large but if the increase is established, 
one can assume that at least a P-wave interaction is present in the 7’=1 state. 
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3) There appears to be a maximum for the bound neutron scattering 
cross-section at about E,= 210 MeV. This could be interpreted as an inter- 
ference between two P-waves. A resonance can be excluded, because no 
maximum has been found for the charge-exchange cross-section in the same 
region. 
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APPENDIX 


Scattering 
No. events No. of prongs | E, AEJE E prongs (*) 0, 

| 

Dublin 
Du 170/11 4 285 63 | 104 133° 
Du 220/11 5 269 51 76 51° 
Du 175/8 1 276 42 30 53° 
Du 177/14 3 295 | 62 56 83° 
Du 175/22 lie 289 49 24 110° 
Du 170/23 3+R 258 84 17 86° 
| Du 194/5 0 292 87 — 725 
) Du 194/15 4 256 78 19 134° 
| Du 185/19 2 252 69 46 95° 
| Du 186/16 0 275 71 — 6° 
| Du 185/23 1+R 238 49 12 42° 
| Du 186/19 2 293 63 76 42° 
| Du 187/4 0 283 84 — 104° 
| Du 185/13 0+R 255 6 — “802 


(*) «E prongs» is the total energy of the stable prongs. 
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Du 188/7 
Du 207/1 
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Du 198/15 


Du 200/16 
Du 200/24 


Du 208/14 


Du 219/5 


Du 189/30 
Du 219/14 


Du 219/2 
Du 219/4 
Du 199/11 
Du 189/15 


| Du 183/9 
| Du 215/6 
©) LE prenes» is total energy of the stable pronss. 
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Scattering (continued). 


No. events No. of prongs E, ABJE, |E prongs (*) 


I 


| @ 
0, 
A = | | | 
Du 201/26 4 274 78 | 125 |. SLA 
Du 198/6 3 275 73 43 82240 | 
Du 191/36 3 271 75 47-50 88. | 
Du 207/13 2 243 66 49. ES 90. 4 
Du 197/14 4 265 86 113 °° | gee 
Du 197/7 2 253 48 72) doom 88° | 
Du 191/32 2 243 77 200 | 58° 
Du 191/34 2 265 65 prese 88° 
Du 210/8 3 283 77 ! 45 |) ARTE 
Du 210/2 2 281 75 179 163 
Du 208/3 LÉ 2h ONE) 2236 68 | 108 | 125 
Du 211/8 3 285 | 72 43: ‘TAG 
Du 214/13 3 284 59 42 | BI de UN 
Du 188/14 0 248 81 — | een 
Du 216/15 1 284 47 104 23° 
Du 218/6 1 236 75 | 9 | 0428 
Du 181/14 4 295 55 40 129° 
Du 202/22 4 265 60 ST CERI 82° 
Du 189/13 1 259 30 | 14500) 30° 
Du 217/10 2 237 71 9 ‘Tab 
Du 217/4 3 266 87 88 V4 
Du 200/12 0+R STE] 7 | «CRUI: 36° 
Du 217/1 1 275 52 | 80 |" * 200908 
Du 208/2 0 290 45 iu 4 18° 
Du 176/42 1 243 56 60 | ieee 
Du 205/26 2+R 260 50 OM 6 | à 
Du 218/2 1 293 89 | SOA 88° 
Du 205/12 3+R Dai 96 ) 17 | 61° 
Du 208/4 1+R 238 60: |) aoe 47° 
Du 200/4 3 | 295 | 66 | 105 D 
Du 203,27 3+R 269 68 dv 49° 
Du 216/16 2 | 300 20 50 ‘ES 
Du 206/16 2 248 64 | 11 | DEI 
Du 219/10 04-2 BON aT | = 59° 
Du 188/23 1 278 35 os o 75° | 
Du 219/7 2 DST CR 58 | 13° > 65° 
Padua 
Pd 7 1 ee a 37 | TR 32° 
Pa 24 0 Mo A RE ee DE rh | 
Pd 58 1+ E Dai pase. ld 9 {7e 
Pd 60 0+R | 210 0 ) ESRI 46° 
Pd 61 Co 242 30 | 70 24° 


(*) «E prongs» is total energy of the stable prongs. 
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Scattering (continued). 


No. events No. of prongs E, |- ABJE | E prongs (*) | 0, 
Pd 73 0+ e ANT 88 = 98° 
Pd 79 4+e 250K. ‘| 85 109 | on 
Pd 80 1 280.0 | 33 47 47° 
Pd 105 sk ls 338 | 78 56 144° 
Pd 114 | 3+R 242. | 81 135 | 52 
Pd 172 2+R 282,0" | 52 CS be 5 90" 
Pd 162 O+R 250 | 36 = 98° 
Pd 209 3 Dani | 64 | 39 | 118° 
Pd 180 1 222 | 50 | = | ZIE 
Pd 228 3 230 | 47 57 | TATE 
Pd 211 0+R Bite |, ‘44 — | 50° 
Pd 248, 2 300 41 10 | 24° 
Pd 250 O+R 200 | 5 = | 54° 
Pd 375 3 260 | "a 23 | 54° 
Pd 396 1 282 87 65 | 21° 
Pd 407 2 235 77 10 | 96° 
Pd 412 DER 220 51 38 | 28° 
Pd 321 MEER 205 | 71 42 56° 
Pd 323 1 ORA 24 | 38 57 
Pd 325 O+R 320,605 | 36 | i 42° 
Pd 264 0 240. | 61 = 74° 
Pd 257 0 255.. | 37 Boe 47° 
Pd 341 0 430 16 | La 34° 
Pd 347 3 st € | 86 | 74 166° 
Pd 444 2 310 | 74 | 115 61° 
Pd 462 BR 23h24 0] — 36 = 
Pd 477 2 325 | 42 30 44° 
Pd 509 3 230... | 52 49 67° 
Pd 512 1 240 | 52 50 102° 
Pd 529 2+R 230 | 52 13 109° 
Pd 531 O+R 235 | 48 = 103° 
Pd 539 4 240 | 69 43 69° 
Pd 590 1 245 | 77 200 120° 
Pd 593 0 200 | 7 =. 41° 
Pd 599 0 Rep SEGGI 48 = 91° 
Pd 615 1 pay kee | 42 = 60° 
Pd 621 5+4R | 430 | 51 82 78° 
Pd 630 1 260 | 20 62 42° 
Pd 643 3 33300. | 88 52 77° 
Pd 652 0 295 | 25 = 45° 
Pd 662 ke, 275 | 40 82 64° 
Pd 664 O+te 2 OSO 48 = 47° 
Pd 669 2 317 73 16 124° 
Pd 694 yey 250 82 59 124° 
Pd 700 4+ R 298 | 54 28 65° 


(*) «2 prongs» is total energy of the stable prongs. 
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Scattering (continued). 


(*) «E prongs» is total energy of the stable prongs. 


No. events No. of prongs E, AE*E | E prongs ©) 0, 
Pd 706 1 lid 308. 58 ae 138° 
Pd 728 1 308 55 66 57° 
Pd 732 2 290 73 131 131° 
Pd 760 3 | 215 62 20 120° 
Pd 774 0+R-+ e | 280 67 — 74° 
Pd 776 - 2 | 250 74 121 95° 
Pd 787 1+ 3B 350 37 13 80° 
Pd 801 2 330 78 29 76° 
Pd 824 1 277 74 | 21 80° 
Pd 844 3 | - 266 65 | 25 80° 
Pd 861 1 | 310 26 | 16 34° 
Pd 864 3 287 71 | 18 101° 
Pd 870 1 244 80 | 2 qTe 
Pd 872 3 230 70 76 92° 
Pd 887 1 317 86 | 135 47° 
Pd 895 0 280 35 | sd, 89° 
Pd 745 1+R 242 60 86 110° 
Pd 725 2 245 56 29 | tox 
Pd 908 1 232 23 7, VAE 
Pd 947 DER 307 89 212 112° 
Pd 966 4 262 52 44 63° 
Pd 970 3 | 212 62 113 32° 
Pd 993 0 242 27 n 28° 
Pd 994 1 210 38 33 66° 
Pd 999 3 285 30 11 43° 
Pd 1017 2 210 55 20 pe ba 
Pd 1031 1 248 27 8 37° 
Pd 1058 1 285 32 31 35° 
Pd 1080 3 230 62 64 97° 
Pd 1120 1 215 51 90 57° 
Pd 1146 0+R fl 242 29 £2 44° 
Pd 1169 2 |. 7228 8 19 31° 
Pd 1172 2 | 400 87 85 122° 
Pd 1193 1 | 340 83 200 79° 
Pd 1204 1 io eg 41 57 27° 
Pd 1208 5 || Bas 65 51 101° 
Pd 1264 PE | 290 36 10 52° 
Pd 1267 3+2R 270 55 48 68° 
Pd 1298 5 305 61 81 107° 
Pd 1308 0+e 260 2 i 33° 
Pd 1314 4+ E 258 62 20 60° 
Pd 1318 1+ E 210 38 3 82° 
Pd 1327 1 292 66 145 90° 
Pd 1335 0+ E 260 44 2” 76° 
Pd 1356 | 2+e 360 79 57 26° 


av | e | 
i | No. events | No. Po prongs | AEJE | E prongs (*) LA 
| | 
| HORS, | 2+R | 77 13 137 
i Pd 1385 | 1 310 es 155 167° 
_Jj Pd 1389 i 1 | 305 33 75 32° 
4 | Pd 1392 | 1 | 260 | 35 24 s2° 
fp nese, I 1 | OT RE = ea 8 93° 
ej ratez | - 2 CD - 7 | 15 121° 
Di Paws | 1 aoe eh ds | 3 91° 
P| Pais. | _o+r i | | — 108° 
» | Pa 1433 CEA 245 15 | 140 108° 
| Pd 1484 1 LENS Re 3 70° 
|__| Pd 1438 4+R suo ‘el - | 57 45° 
ie - 1 | 270 | 63 .| 12 70° 
| PA 1474 Nic sala "6 | sE SES | 52° 
3 i (*) «E prongs» is total energy of the stable prongs. 
Charge-Exchange 
È 4 3 | Enegy of | 8 of 
— a St pense si; fastest protons | fastest protons | 
Dublin 
Du 220/8 e 2+ RR | 182 1 = 
Du 177/3 | 3 | 282 52 62° 
Du 175/11 2 +R 260 107 38° 
Du 193/3 4 293 12 63° 
Du 175/2 | 2+e 282 4 107° 
Du 185/24 4 251 29 46° 
Du 186/21, 3+ = | 285 45 29° 
Du 186/22 | 1 | 290 5 149° 
Du 187/11 3 | 262 38 60° 
Du 191/5 | 2+R | 278 32 78° 
Du 193/2 | 3: 234 15 53° 
Du 202/6 5 239 96 29° 
265 | 8 138° 
# 252 | 9 = 
ee S 86° 
| 239 | 9 108° 
| 272 | 165 63 
| 244 | 74 | 93° 
| 263 47 7 
269 95 49° 
263 | 7 | Tr 
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Charge-Exchange (continued). 


\ 


HAE 


GRILLI, L. GUERRIERO and G. A. SALANDIN | 


der lab] 


j—~Energy of 0 of 

No. events No. of prongs LA fastest protons | fastest protons 
Du 198/4 6 281 15 48° 
Du 205/1 2+ E 275 91 52° 
Du 198/9 3+ E 290 80 651° 
Du 204/5 1+ Rk 253 3 104° 
Du 202/4 2+ Rk 287 53 30° 
- Du 219/11 5 280 51 108° 
Du 207/21 3 281 27 50° 
Du 191/27 3 271 16 ‘1139 
Du 219/3 3+ fh 287 112 Sie 
Du 205/10 2+ Rk 239 116 250 
Du 218/1 il 270 10 147° 
Du 201/35 1 285 86 203 
Du 183/24 1 236 115 32° 
Du 184/31 Il 289 93 30° 
Du 213/9 1 242 92 20° 
Du 176/38 2 243 25 20° 
Du 216/17 4 250 23 63° 
Du 177/20 1+ Rk 250 82 59° 
Du 200/35 3+ E 295 46 08 

Padua 

Pd 18 3+R 255 54 30° 
Pd 28 6 270 36 38° 
Pd 46 2 277 45 48° 
Pd 54 3 210 12 — 
Pd 55 4+ Rk 275 Hl — 
Pd 66 0+e 350 — — 
Pd 84 2 307 170 57° 
Pd 88 l+e 242 150 Sue 
Pd 103 1+ 8 262 12 _- 
Pd 126 2 225 58 025 
Pd 210 3 262 199 10° 
Pd 203 2+e 250 | 160 16° 
POS 2 238 80 69° 
Pd 451 3+R 260 168 24° 
RATA: 2 310 50 se 
Pd 510 5 440 127 24° 
Pd 519 5 330 165 144° 
Torok Mayra 3 300 80 67° 
Pd 432 2 282 245 42° 
Pd 557 2 200 72 69° 
Pd 632 1+£ 318 180 98° 
Pd 709 5 310 27 88° 
Pd 778 5 300 40 54° 
Pd 791 1 240 16 753 
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Charge-Exchange (continued). 


| No. events No. of prongs | E; er oi | N aces 
fastest protons | fastest protons 
Pd 823 2 270 240 55° 
Pd 853 SAR 245 13 34° 
Pd 555 CdA 262 70 | ire 
Pd 755 2e 210 8 — 
Pd 340 1 | 280 225 | 15° 
Pd 907 342 OR 286 35 | 120° 
Pd 913 4 210 34 3° 
Pd 980 1e 245 10 02° 
Pd 1009 4 262 19 2 
Pd 1014 SR 300 | 45 | 43° 
Pd 1068 5 210 52 | 14° 
Pd 1084 3 330 11 | = 
Pd 1102 1 250 TO. | 49° 
Pd 1103 5 272 50 200 
| Pd 1143 SEL 325 19 125° 
Pd 1144 1 270 | 130 43° 
Pd 1159 3 275 | 29 177° 
Pd 1170 IERI 215 | 80 35° 
Pd 1199 Le 380 133 39° 
Pd 1224 1+R 282 230 23° 
Pd 1280 2-7 278 54 92° 
Pd 1281 | 5 260 56 86° 
Pd 1309 i; 300 97 103° 
Pd 1330 2 270 43 67° 
Pd 1357 4 262 74 35° 
Pd 11360 | 2 | 280 19 45° 
Pd 1480 | 1 302 27 99° 
Pd 1527 2+R 350 18 36° 
RIASSUNTO 


Sono state analizzate le interazioni relative a 185 m di traccia di mesoni K° osser- 
vati in emulsione nucleare. Si sono ottenute informazioni sulle reazioni di scattering 
K*-p e K*-n e di cambio carica. Basandosi su tutti i dati in emulsione raccolti nei 
vari laboratori, si è determinata la dipendenza dell’energia delle 3 sezioni d’urto fino 
a 350 MeV. L’analisi in fasi mostra che verso i 300 MeV il contributo dello stato T—0 
è dello stesso ordine di grandezza di quello dello stato T=1. All’interferenza tra onde P 
di questi due stati si può attribuire aumento della sezione d’urto di cambio carica 
e la diminuzione della sezione d’urto di scattering su neutrone. 
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Partial Wave Analysis of the Production of Boson Pairs. 


S. CIULLI (*) and J. FISCHER (**) 


Joint Institute for Nuclear Research - Dubna 


(ricevuto l’11 Gennaio 1959) 


Summary. — Partial wave analysis for boson pair production on nucleons 
is made. The corresponding angular operators, which characterize the 
spin and angular dependence of the S-matrix, are expressed with the help 
of the Legendre polynomials and tabulated. Simultaneously a general 
method of calculating the angular operators is given for processes con- 
taining more than four particles, in which cases the straightforward 
method leads to very lengthy calculations. 


4. — Introduction. 


Boson production processes as pion and photon pair production, radiative 
pion scattering etc., have been intensively investigated during the last years. 
The theoretical bibliography on this subject may be divided into two parts: 
one containing pure-theoretical treatises, which solve the problem by means 
of the Chew-Low method or dispersion relations, and the other consisting of 
works which use some semiempirical postulates concerning the dynamical 
character of the process, such as assumptions on the existence of resonant states. 

For both these methods, it is useful to study the angular and charge de- 
pendence of the S-matrix. In the first case, one obtains in this way a set of 
independent equations for the energy-dependent coefficients; in the latter case, 
one can write directly the amplitude of the assumed resonant state as a function 
of angular and isotopic variables. 


(") On leave of absence from the Institute for Atomic Physics, Bucarest. 
(”) On leave of absence from the Physical Institute of the Czechoslovak Academy 
of Sciences, Prague. 
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To perform this angular and isotopic analysis, it is sufficient to know the 
general laws of conservation, without assuming anything concerning the dyna- 
mics of the reaction. The result of such an analysis is a set of orthonormal 
polynomials, which can also be used for the phase analysis of the experimental 
data. This is especially important for processes with an interaction character 
which is not very well known. 

In the present paper, the ordinary space structure of the S-matrix is studied 
for processes which may be deseribed with the help of the following formula 


(1) b+N—)b4N, 
t=1 

where N denotes the nucleon and b, b; bosons of any kind. For processes of 
the type b+N—b’+N this analysis has been made by RITUS (1). However, 
the straightforward generalization of his method to the case of more particles 
leads to lengthy and cumbersome calculations. For instance, for n= 2 in (1) 
almost two hundred terms must be calculated and for higher values of n the 
calculations are practically impossible. We give a method for simplifying this 
procedure, namely by reducing the production of n bosons to that of n—1 
bosons. We demonstrate it for n= 2, but it has general validity (n= 2 is 
chosen only to deal with simple formulae). For the same value of n we have 
calculated the explicit form of the angular operators. The resulting polyno- 
mials are given in the Tables I, II and III for different coupling schemes of 
the final angular momenta. 

Readers interested only in the practical use of these Tables can omit 
Sections 3 and 4, which deal with the calculation formalism of the angular 
operators in the boson-boson and boson-nucleon coupling schemes respectively. 


2. — The angular operators. 


Let the initial state of our system be characterized by the total angular 
momentum J, its <-component M and by the eigenvalues (i) of a system of 
quantities which together with J and M form a complete set of commuting 
observables. Then, the initial state will be described by the following ket- 
vector 

J M(i)>. 


Similarly, the final state will be described by 


TMi); 


(1) V. I. Rirus: Zu. Eksp. Teor. Fiz. SSSR, 32, 1536 (1957); English in Soviet 
Physies (JETP), 5, 1249 (1957). 
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where the symbols have an analogical meaning to that of J, M, (i) respectively. 
The operator which describes the transition |JM(i)> —|J'M'(f)> is defined by  « 


= 
x 


FMP TM . 
Consequently, the S-matrix can be expressed as follows 
(2) S= Ya(J'M'(f); IM) |I Mf) TMi) |, 


where a(J'M'(f); JM(j)) determines the « breadth » of the corresponding chan- 
nel. Doing the partial wave analysis, we shall be concerned only with the 
angular momentum part of the channel labels (i), (f) and we shall omit any 
explicit reference to other quantities such as energy, isotopic spin etc., which 
would be included in the a-coefficients. Due to the three-dimensional rota- 
tional invariance of the S-matrix, the right-hand side of (2) may be written 
in the following form: 


aI) STM M<TMA |, 
where ; 


(3) SJ MT M) | = FINO) 


M=-J 


are the so-called angular operators of the reaction (see (1)); they determine 
the angular and spin dependence of the S-matrix. They are orthogonal with 
one another and are normalized in the following way: 


tr FHI) FIMO) = (CI + 1) Our dd 


over 
spins over all 
angles 


if the initial and the final state eigenfunctions are normalized to 1. 

Before calculating the explicit form of (3) for reactions belonging to the 
type (1), let us remark that if we know the angular operators for b, b; (see (1)) 
having spin zero, then those for spin one may be obtained by operating with: 


il 
fos for the electric multipole , 


V Ix (Ip +1) OF 


— i © ; 
-—— |k | for the magnetic multipole and 
VI (Ii il ok 


k for the «longitudinal » multipole , 
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on the eorresponding angular operators. Therefore we calculate only the angu- 
lar operators for b and b; having spin zero, and we do not suppose any defi- 
nite parity for them. (So, if some of the b, b; are 7-mesons; the corresponding 
angular operators are obtained by postulating the internal oddness of these 
particles.) In this particular case, the set (i) of the initial eigenvalues consists 
of the spin o — + of the nucleon and the orbital momentum / of the ingoing 
boson. The symbol (}), in the two-boson case, represents by itself either 


(1) o = } ’ L,, l, ’ L where È = Ì, + Ì, 
or 
(II) I) rer TRAME À More “A Se 4 i, i 


Both these coupling schemes must give the same information about the pro- 
cess, if transitions from all possible initial to all possible final states are con- 
sidered, 4.6. if no one from the a-coefficients is neglected. Mathematically, 
this is an evident fact because of the completeness of the orthonormal set of 
the angular operators both in Case (I) and in Case (II). In practical calcu- 
lations, however, we deal always only with certain non-complete subsets of 
these sets, having to restrict ourselves to a finite number of terms in (2). 
Then, the two coupling schemes are no more equivalent and we must choose 
that of them which gives a better approximation to the ideal case. This de- 
pends, naturally, on the distribution of the a-coefficients in (2): if only few 
of them are large and the others small, then it is sufficient to consider only 
a small number of channels fer obtaining a true picture of the reality. If, 
on the contrary, we deal with such a coupling in which all or many a’s are 
large, the calculations approach very slowly to the real case. Thus, this 
practical reason gives a criterion in choosing the coupling. 

Besides this practical argument, from the physical point of view, the fact 
that in a given coupling scheme only a small number of channels are important, 
gives a deeper insight into the understanding of the nature of the process. 
Indeed, the existence of such resonant states cannot be predicted by means 
of general group-theoretical methods, but is in a direct connection with the 
dynamical properties of the interaction. For instance, it is to be expected 
that the isobaric state (3, 3) of the nucleon will play an important role in 
processes in which at least one 7-meson is present, as it is the case for nu- 
cleon-nucleon collisions and for elastic scattering of pions (2). 


(2) We refer the reader to the clear discussion on this topic contained in the paper 
of R. F. Pererts: Phys. Rev., 111, 1373 (1958). 
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In general, we can say that the choice of the coupling scheme is determined 
by the interaction: if in a system of three particles, a certain pair interacts — 
much stronger than other pairs, it may be considered approximately as a 
shortly-living subsystem with spin-equal to the vectorial sum (according to 
the vector model) of angular momenta of these particles. Naturally, by spe- 
cifying the spin eigenvalue of this subsystem we determine, simultaneously, 
the coupling scheme of the angular momenta. This situation is quite analogous 
to that in an atom, in which the angular momenta may be coupled either 
according to the (ZS)-scheme or according to the (jj)-scheme. 


8. — The boson-boson coupling. 


In the case of the boson-boson coupling, the final state may be written 
as follows 


z 
|JMUILI) = > CMH gah |} 


LM—p 2 
p'=-% 


where 


Qu - (q, 1) de Ch VELE 2H — n° FUSE wD, LE PAT 


MES 


is the orbital eigenfunction of the outgoing bosons and |4u'> is the spin eigen- 
function of the nucleon. If we write the initial state also as a Clebsch-Gordan 
combination of the angular momenta of the ingoing particles, we find that 
the angular operator (3) has the following form 


(4) | FULL) = > 3 Ly cha Gul, 
pet p=-h 

where 

(5) Lay chi #2 Ca ia Or par ER #(Dr rx 1M- Es PS 


will be called the orbital operator and 


fi i! 
for p= È == 

(6) di Bi i 
| de, tor w= pp SE 


PARTIAL WAVE ANALYSIS OF THE PRODUCTION OF BOSON PAIRS 269 


will be called the spin operator (¢,, 0,, o. are the Pauli matrices); p, q and r 
are unit vectors parallel to the momenta of the b, b, and b,-particles respectively. 

Comparing this with formula (5.II) of (!) we see that the angular ope- 
rators (4) have the same form as those for the elastic scattering of (parity- 
less) pions on nucleons; the only difference is that the orbital eigenfunction 
of the final pion is replaced here by the Y-function. Therefore, all calculations 
are analogous to those for the elastic scattering, namely, the sum over 
M=-_-J,...+J in (5) can be reduced, by choosing the z-axis parallel to p, 
to one term: 


L "+ Ick 
Ce Er ICRA S Ba u (q, r) re a) 
J 


=] 
— 


(because M — u — 0 in this co-ordinate frame). By inserting this expression 
in (4) we obtain the form of ¥ in a special frame and, using the rotational 


invariance of 7, we can write it in an invariant form 


Let us perform the calculations in some more detail. The first Clebsch-Gordan 
coefficient in (7) is different from zero only for J= L + 1 and the second one only 
for Ja +43. Un u= ++ and p’= ++ so that there are sixteen values 
of 2, for given 2, L, ly and 4. Let us remark now that the magnetic index py — u' 
in (7) has only the values 1, 0 and — 1, so that the @-function may be expressed with 
the help of the operator 


in the following way 


(8) VI L+ 1) @ ir = = LV 


where Lc TA + De Now, according to (6) and (7), the four terms on the right- 
hand side of (4) may be written as follows (we take, for example, the case J=L+3= 


=1+)): 


L +1 1 + 0, 
DE DEEE mls (gr) for u'=u= +4, 
/ 20 +1 2 
(9) and 
1 = O, + 10» 
ee fg rh = ; for w=—u— +} 
ASE A 2 
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They depend, evidently, on the co-ordinate frame chosen. In order to write them 
invariantly, let us define the vectors i 


P,=r-ppr,  P,> (pry, 


which together with p form a carthesian system and are, in the p|e-frame, parallel to 
the x, y, z-axis respectively (the z-axis is chosen in the plane of the vectors p, r). Their 
lengths are |P,|=|P,|= V1—(p-r)|?, |p|= 1. With the help of them, we can write 
(9) in the following invariant form: 


1 


|2P,|}V2L+1 


where P, = P, + iP,. Taking into account that 


P--LYVnP,:0, 


PL P.,-c+P,-LP_-c ; a 
2|P,| SLA 


and that Poe 0 we obtain the following form for the angular operator: 


(10) F=— 


In this form, the angular operator is still independent of the number of outgoing par- 
ticles (for more particles it is only necessary to write more arguments behind @);;). 
For the two particle case, we can write 


li 
TAU, fh, eed 
(11) iy Cre 4 NV ey ieee 
dati 


For simplicity, we restrict ourselves to 2,= 0 and 1, J, being arbitrary. We obtain 
for — 0 


eae 4/21, +1 
VA de de 


Ys = Pi, (pr); 


and for 1,= 1 


AIRONE V3 
DIO 


n A, + 1 
11, = . 2 À I È 
1, feu LU + De [gr]D.(p.r). 


via, | 
HE Teves (—lLp-qP,(p:r) —[pq):[pr1Di,(p'r)) ; 


(+ Dp-qDi,(p-r) —[pqlprlD;(p-r)) » 


where D,,(p:r) is the Legendre function of the order /,. 
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By inserting them in (10) and in analogical expressions for J=L43=13 3 


J=L—-3=1-4 we obtain sixteen different angular operators as given in Tables I 
and II. 


4. — The boson-nucleon coupling. 


In the case of the boson-nucleon coupling, the angular operator is given 
by a formula which is quite analogical to (4): 


+ 4 


(12) FIhbj}) = > È Futa su')Gu|. 


= 3 


However, we can easily see that this form of writing, in which the orbital 
and spin quantities are separated, is very disadvantageous in this case. Indeed, 
£,, in (12) have the following form: 


u'u 


J 


(SAME > SIL (lp: 
M=-J m=-j 

If comparing it with (5) we see that (13) contains two summations (namely 
over M and over m), while (5) has only one, namely that over M. The second 
Summation in (5) is involved implicitly, through Y, which must be expressed 
through Y(q) and Y(r) according to (11). For this reason, even if the sum 
over M is eliminated in (13) by suitable choosing of the co-ordinate frame 
(as it was done in (5)), the sum over m remains. This circumstance, naturally, 
complicates considerably the calculations. Therefore it is more convenient to 
write (11) in another form: 


4 


(14) F(T lj $14) = > > Fray 13M E ul, 
m=tTÎ u= + 
where 
(15) de 2M — LEO LA ve lf r). Te) 
M=-J 
and 


\jml,3 3) = a 5 OS nee Y Y,, mp (q) | 


p'=-$ 


It is clear that all the four summs will subsist in (14), namely the two 
« small » (over w and y’) (*) and the two «great» ones (over m and M). The 


() One of them (u’) is contained implicitly through |jml,$ 
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last two are responsible for the existence of a great amount of terms, which 
make all direct computations impossible. Nevertheless, (14) has the advantage 
of an explicit summation only over m, M being included in &,,,, (see (15)). 

To get rid of the summation over m, we have to take the vector q parallel 
to the z-axis; on the other hand, the expression for #,,, reduces to one term 
only in the co-ordinate frame in which p is parallel to the z-axis. Since #,,, 
is not invariant we must, having found its form in the plz-frame, transform 
into the g|z-frame, in which (14) may easily be calculated. For this reason, 
we shall investigate the transformation properties of #,,,. Let us consider 
the following angular operator 


R (J1,j13) ni 2% my lI™>< Sp | - 
3 


The only difference between (16) and (14) is that |jm> in (16) is a pure spin 
function corresponding to spin j, while |jml}) in (14) is a combination of a 
spherical harmonic Y,,,_,,(g) and a one-half spin function |3u'). The ope- 
_ rators 


dm p| 


describe the conversion of the nucleon into a particle with spin j. Thus, the 
R-operator defined by (16) describes the elastic scattering on a fermion the 
spin of which changes from } to j. 

The method consists now in the following: from (15), we find the form of 
Kay in the P iz-frame and inserting it in (16) we obtain the explicit form 
of R. Since Ris invariant, it has the same form in the q|z-system. The form 
of X,,, in this system may now be obtained by means of a trace operation 


on the product of R with the spin transition operators ((jm>)y|)}, written 
in the gz-system: 


| R = YA) limp» Gutelp)| 


may 


(17) A 
| 2,619 = te (F-uela)Gmela)1), : 


where |}u(z\q)» and <jm(z\q)| are the spin functions in the q|z-frame, for 
spin 4 and j respectively. In order to write easily these transition operators 
in different co-ordinate frames it is convenient to express them as linear com- 
binations of the spin-tensors of Racah, which, as it is known, provide a basis 


for the corresponding matrix algebra. Moreover, they have the advantage | 


| 
| 
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that their transformation properties are known, so that they can easily be 
written in the desired co-ordinate frame. The spin-tensors are defined as 
follows: 


(18) TE (j, 4) = (IAC 5 E— | 


and transform according to the irreducible representation of weight 2/+1 of 
a 
the three-dimensional rotation group (in our case [=j+4 or I=j— 3. 
Let us perform this calculation for j=} and j= 3, for which values the 


Clebsch-Gordan coefficients are tabulated (see (*)). We obtain, for j=4 
Hey ak 0, + 10, 1 
E È st ae ade 0(1 1) — 
(194) PES) == (5 1} MAL, 1($4) = + 9 ’ T2 (82) =F Oz; 
and for j= 3 
MEO Gx 0 05-46 
A LIBRO. A V2 ..0 a 11 0000 
ITR re TONS Sie es a ia ee È 
Ana oat PA a |’ ADIEU DE 
OO ON <6 0 v3 
0 1 1 0 0 0 
st 0 0 3 ig Ea ees i, LIV2 0 
19b) nas Pe a 20 
Oe i DAR 0 po CA I 
OW) 0 0 \ 0 0 
OO CORO 
ae x | ANT) 0 k 0 0 
om noi ie regi L + 
2004/80] 0 0 
DRE 49; 


In order to be able to write them in different co-ordinate frames, it is con- 
venient to make use of the carthesian components 7;, 7,., rather than the 


(8) E. U. Conpox and G. H. SHortLEY: The Theory of Atomic Spectra (London, 
1935). 
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circularly polarized ones. They are connected by the following relations: 


aa de 
(20a) perzq VE? Ta 


(208) EE INIZI 


Ria 1 
122 REV (2T,, — Tra — ee 


43 Bo 


(see Appendix). Now, if we use the relation inverse to (18) we can express, 
with the help of (20), the spin transition operators |jm><4| in their tensorial 
form. We obtain (omitting the labels 7 and 3) 


for 7=4 
| lepoai= Lovin, 
V2 
(21a) 1 
EDF = et). 
and for j — à 
PERV CA Pl ae seer Ea (Pees 
aL 1 
1 i|/— = = DLE, Je o) aa 
His So ) 4/2 3 
(21b) 
V3 RE Tu 
Fb) = 5 Maud 4%) 2V2 (T, il); 
FEAR =F 37 To) + Ton 


Now, we must compute Æ,, both for j=3 and j=}. For 4 =#%, there are 
eight possibilities of combining 1, with j and / with o=4, namely: 


TERRE VE, 


T=h+3=14}; 
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For j= +, there are the following four possibilities: 


Let us perform the calculations for an example, say J=1,+3=1+ 4, i.e. l=1,+1. 
In the co-ordinate frame in which e|p and the x-axis is chosen in the plane of the 
vectors p, r, we have, using the Clebsch-Gordan coefficients: 


1) + 2) es —_ 
2+ 1 — 
n Ta Loan Vas ns MVA+ 


This expression is to be multiplied by 


1 1 8 1 
Ta ch dI) ae sl 2 (el iT,) 


A 


Since, however, 


È h+2. 


Re ri Di (p:r) 
ae 4m 4/ (21, + 2 2)(2ly + 3) ; 


is to be multiplied by 3(7,, — iT,.) —3$V3(T,—iT,) and since r.=r,=r, in our 
system, the sum of these two terms gives 


i> 
1 1, + 2 F 83 
a Di lpr) rd (iz. +|/£ r.) 
470 4/20, + 2)(21,+3) # 


This expression can be written in the following invariant form 


1 1, + 2 1 (E CE | 
lea pira. 
er fe ere 7 21/2 È 


In a similar way, terms containing Ayg=}, 1343: Au3,33 are calculated. Sum- 
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A 


ming all these contributions, one obtains the following expression for A: 


CS 1/(4x) 
R =. ($7) —— {—V6( + 2), +) p-P Rip) + 


V/ (2l, + 221,4 3) 


(GP iP, -T-p) Di(p:r)+ e 


Similar expressions were obtained also for the other seven cases of j= 8. 
TS 
By the same method, also the A-operators for j= 4 were obtained. 


To find the 7-operators, we make use of (17). It is easy to see that the 
orbital operators £2,, may be written as follows 


Fan = Tl ca CM VER FI te Cu) Ge |). 


The trace operations can quickly be performed using the formulae 
OT) = 0a 
tr (TE Tu) = $(Oix0j2 + a O5) 
MS eat Ore 


These relations have only a restricted domain of validity, but for our purpose 
one can use them without supplementary precautions (see Appendix). 


For instance, for 7 — 3% (4 =1), J=1,+3=1+4 we obtain the following result, 
if we are working in the q||z-frame (*): 


{ (le + 1)(le+ 2)p-qDi, + 


+ (l,+ 2)ip-q P,-q —P.-q)Pi,— iq-P,q:P.Pi}. 


The other orbital operators have a similar form or are somewhat more complicated. 
By multiplying £,,, with |} u'><$m| and summing over y’ and y, we obtain the cor- 
responding angular operator. The resulting angular operators are tabulated in Table I 
and Table III. 


() The result is dependent on the co-ordinate frame. 
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All practical calculations of Section 3 and Section 4 were checked. More- 
over, the normalization and orthogonality of each angular operator was ve- 
rified by direct integration. i 


5. — Conclusion. 


We must mention that the calculation of the angular operators is only 
the first step in investigating the S-matrix of a given process. To do the se- 
cond one means to perform the same analysis in the isotopic space. The method 
of obtaining the corresponding « isotopic angular operators » is quite analogous, 
the only difference being that if one or more photons are present, the S-matrix 
is no more a scalar but rather a sum of a scalar and the third component of 
a vector. This means that instead of each angular operator in the ordinary 
space, there are one scalar and three vectorial operators in the isotopic space. 
So, the number of operators is greater, but the calculations are not more com- 
plicated because the whole procedure is quite automatical. 

The third step consists in constructing the theory of interaction. This 
problem must be solved by other means than the first two, for instance, with 
the help of dispersion relations or of the Chew-Low equation. However, it 
should be mentioned that for the complete analysis of the reaction the first 
two steps are necessary, because only after having done them, one can com- 
pare the theory with the experiment. Indeed, if the angular analysis gives 
the explicit form of the angular operator HI, (f), (i)) and the form of the 
S-matrix is obtained by use of the interaction theory, the experimentally 
observable coefficient a(J(f)(i)) is defined by the theoretical formula 


a(J, (©) = te [FI © GS, 


where tr Î means integration over all continuous and sum over all discrete 
variables. In the isotopic space, the situation is quite analogous. 


We would like to express our gratitude to Professor NING Hu, Professor 
S. TireicA and Professor YA. A. SMORODINSKY for many illuminating discus- 
sions. We thank also Mr. V. G. Erimov for helpful conversation during this 
work. 


18 - Il Nuovo Cimento. 
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Tables of the angular operators. 


These Tables contain angular operators for processes of the type b+ NV > N+6b,+)b, 
(where b,, b, and b are zero-spin bosons without definite parity) for 1,=0, 1 and 2, /, 
being arbitrarily large. 1, 1, and 4, are orbital quantum numbers of b, b, and bd, re- 
spectively. L is the intermediate quantum number in the case of boson-boson coupling 


in the final state, (È = i+ i), j is the intermediate quantum number in the case of 
boson-nucleon coupling in the final state, ( j= Pie) P. q; r are unit vectors in 
the directions of the momenta of b, b,, b, respectively. For practical reasons, we make 
use of an orthogonal set of (not normalized) vectors P,= r —p pr, P,= [pr]. p. 

Pi, denotes the Legendre function depending on p:r=cos 0; Di, and Di, are 
derivatives with respect to (p-r). 

If one of the bosons is a vectorial particle, the corresponding angular operators 
can be derived from ours by means of the following operations: 


1 
n for the electric multipole, 


Vi, (1,4 1) OF 


for the magnetic multipole, 


— d 
= es we le pe 
Vk&+1| Ok 

k for the longitudinal multipole. 


Here k is the unitary momentum of the corresponding boson. 


TABLE I. 
(I=L+3=143 4=0 Len 124 
|J=b+3=1+3 4=0 fat I=) 


= 


F = (dr)? {(l, + DD, + io-P, Pi} 


ie TE 
ba ei 


F = (40) {— (+ Do-pDi, +0-P, Di} 


— 
N 
Pe 
s bi 
a 
(SEI 
I 
ti 
we eH 

I 


x 
= 
Si; 


[I=L-3:=1+4 k=0 Deh d=k_-1 
|J=4-4=144 4=0 .j=3- 1=h-1 
F = (42)? {lo -pQ,—o-P,D,,} 
(I=L-3=1-3 4=0 L=k I=), 
[\q¥=4-$=1-4 4-0 G=3 I=] 


L 
ID 
F = (4n)-#{l, Di, —io-P, DI) 
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Le] 


DD PP TEE L=1 L=Lh+1 I=1+1 


1/8 


F = EG, 4. 1)((, + 2)p-q + io-[pg) Di — 


V (lp + 1)(24 + 3) 
— (1, + 2)(P,-q—ie-P,p-q) —io-[pq]p-r) Di, —io-P,P,-q Di, 


=_= + INC + 3)o-pp.q -0c-q)D, + 
Vile + DE +3) Me MT ae 


+ (@+ 2)(6-pP,-q + oP, p-q) + (6-q —o-pp-q)p-r) Pi, —o-P, P,-q Di, 


Je tee =) DSRt1 124, 


ott A] 


( + 1)(L0 : ‘a +6:q) Di, + 
PRIT CMS 


+(@+1)c-pP,:q — (l,+2)6-P,p-q —(6-q—o pp ‘q)p:r)P:,+0-P, P,-qDi, 


SRI Ft) Lil 1141 


F = (4)! 


V3 
QUEI) (2,, + 3) 


{( + 1)((% 1)p-q —io‘[pq)) D,, — 


— (+ DP, -q + il, + 2)0-P, p-q + io-[pqlp'r) Di, + io -P, P,-q Di} 


ete LL b= 121, 


1/3 


7 = (42) ——— {(i(l, +1) P,-q—o-rp-:q+o-qp:r) Di, —o-P, P,-qDi,} 


~3 { 
NC 


LL LL. 1—1+1 


3 = (4a 8 


n 


{((+ 2)c-p P,-q+6-[qr]) Di, —o-P,P,-qDi} 


MER) 
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44 , My 4 a re r 1} 
ò e > » + f Pad avin he A 
pet LES put. at EN eo, + i i ù AI £ 


fes 280 LR ONCE be PI ent and 3. riscuee 


Taste II (continued). \ 


4 ILlI<4=144-L=4 fea V= LA 
È. A — 3 
9 = (4094 7 {((L—1o-pP,q_-0- ae P,+0-P,P,-qPi} 
ie V/1,(l,+ 1) 
os 
ae: CIO = AT it 
“1110 2 | 
a Lx 7 = (42) Mr ar {i,P,q+orp-q—o-qp'r)D:,+o'P,P,-qPi} 
27 22 ; 


Sa T44wÌ43 al Dal Ali hi 


= (ul 


Fie the 14 La bel Ve 1] 


4 MICI 
7 VL(2L — 1) 


+ (Lo:pP,;q (I, 1)c-P,p:q+(0:q-0-pp'q)p'r) Pi, -0:P,P,q Pi} med DL 


{h((l,+ l)o- pp-q —o:q) Di, + 


-J=L-3=1+443 h=1 L=4-1 h=h-2 


V3 
= (Artt a L,—2)o-pp- 7 
ÿ = an) vira Ta ((l,— 2)o-pp-q + 0°q) Pi, + 


Sa Ter fe ne Di, +0-P.P.qDi} | 


Ione LD on Vin fn LT ee an LT 


3 = (401 Sa {1 (hp: q—~ta-[pql) Pi, — 
22 


— ((l,—A(P,-q — is P,p:q)+ io lpglp:r)D',+io-P,P,- ap} 
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TABLE III. 


F=- (an) KE + lo-q Pi, + i(P,-q —io-[P.q]) Di} 


Fate ts ye Se t= at I 


= 


F = (4x) ={(& + 1)(p-q — io-[pq) Pi, — (P.-q — io-[P,ql)D:.} 


eee A ae] 


F = (4n)-1{l.( p-q —io-[pq]) Pr, + (P.-q — io-[P.q)) Di} | ig 


Pee tt bal = t=] 


F = (42) {Lo qDi,+iP,q — io [P,q]) Di} 


“ 

7 = An ——__"___— {(I, + 2)(,+ 1)(2p-q+ io-lpq])) Pt È 
OV (ly + DCR + 3)” PI Payer, à 

È 

+ (b+ 2)(i{o:P, p-q+o-p P,-q) — (2P.-q + ic-[P.g))) Pr, - 

Te i(o-P, IT o-P,P,-q)D:} | ; : 

i 

MERE = ges 11,42 ; 
Re = {= (+ 2)@ + DBo-pp-gq —F-q) Pi, + à 
2V (lp +1)(2l.+3) ; x | à 

+ 2% + 2)(6-P,p-q +o-p P,-q) Di, —(6-P.P..q —o-P,P,-q) Di} | sa 


DERE LEE L=l j=s 1=L 


2v/1,(2l, + 3) © DETENTE 


—(3(6-P, p-q + o-pP,-q) —i(2,+ 3)(2P,-q + io-[P,gq))) Pi, + 
“| 3(0-P.P.-q —e:P,P,:q)Pi} 


3 = (42) 
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[is ee 1 ee ee LEI 
À 1 
= (4)! —— {l(l,+ 1)(2p-q + ic-[pq] Pa 
F ai) PF PQ) Pi 
— (3 + 2)(6-P,p-q+o-pP,-q) + },(2P.-q + io [P.q])) Di, + 
Sa 3i(o-P,P,-q sla o-P,P,-q) Di} | 
Eyes ee CT LD a eg 
= 1 
= (47) =? —___{_1,,+ 1)(2p-q + io-[pq]) Pi, — 
3 x nh PY Pq) Pi 
— (3i(1, — 1)(6-p)P,-q +o-P,p-q)+(+ 1)(2P,-q + io-[P.q])) Di, — 
— 3i(o-P,.P,-q+o-P,P.-q) D;} 
Sp io ES e LE Gay fee ener al 


-{L( + 1)(36-pp-q —o-q) Di, + 


+ (3(0-P,p-q+o-pP,-q) +i(2 — 1)(2P,-q +io-[P,ql)) Di,— 


— 3(6-P,P,,-q —o-P,P,:q) Pi} 


F=l-3=144 (L=1 j=$ h=h=-2 
à | 13 
F = (4n)-t ————___{, (1, — 1)30-pp-q —o-q) Di, + 
Naga E PP'q q) Pi, 
+ 2(l,—-1)(6-P,p-q + o-pP,-q)Di, + (6-P,P,-q —o-P, P,-q) Di} 
To=35 4, Le lege? i t=h—1 
Ila {- L(l,— 1)(2p-q + io-[pql) Pi, + 
= a PI . 10° 
2 CL — 1) i 1) 202 70 q PW le 


+ (&—1)(i(o-P,p-q + o-pP,-q) — (2P,-q + io [P,q])) Di, + 


+ i(c-P.,P,:q+0:P,P.-q)Pi} 
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APPENDIX 


Relation (18) is a unitary transformation connecting the spin transition 
operators |jm)<3u| with the circularly polarized spin-tensors 7’+? and 7°. 
We see from it that the circularly polarized spin-tensors are orthonormal: 


(A.1) (CTP) TS) = di ,dee . 


Further, we can use (18) for recognizing the transformation character of the 
spin transition operators |3m><4u|. For this purpose, we express 7? and T1 
through the corresponding carthesian components 7,. and 7,. The relation 
between 7! and 7, is well-known: 


e Lar 
(A.2) À V2 ( at 1) ’ 
| 719 — T, 


To find the relation between 7? and 7,,, we express T? as a Clebsch-Gordan 
combination of two independent vectors. We obtain 


| 1 
T2=? = 3 Ud, 


Con i i 
(A.3) = 3 (uv, + Wz) — 5 (U,0-+ uv) , 
1 1 2 
peso UE UV a U,0, + IL UD, 
and 
i 
sa (UV + U,0, + Ue) - 


This is a unitary transformation between the circularly polarized tensors and 
the symmetrized pairs. The corresponding carthesian symmetrical tensor is 
defined by 


f ISAIA 


AA 
oe | Toy = Typ = Di (Wy + U,V2) à + 
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To make the transformation (A.4) unitary, we define also 


(A.5): 


In this way, 7,, will be expressed through 7? and 7°°= T by a product U 
of two unitary transformations. However, being not contained in the rela- 
tions (18) for 7 = 3, 7 cannot be expressed through the spin transition oper- 
ators |3m)<3u|. Therefore, it must be constructed artificially, with the help 
of a supplementary vector, p say. For instance, let us take 7 =p-T. Then 


(A.6) tr (T+T?*) = 0, tr(T+T)=1. 


Hence, according to (A.1) and owing to the unitarity of U, T will be an 
orthonormal set of matrices 


[aeree 
(A.7) 1.6. 
| tr (TT) = #(0;%0 5 + O10 jx) « 


Combining now the three transformations (A.3), (A.4) and (A.5) we obtain 
the explicit form of the matrix 


(wx) (yy) (22) (ay) (y2). (22) 


1 | i 
e 0 7 0 0 
i 1 
0 0 0 GET E 
V2 V2 
T 1 2 
+ = 0 0 0 
V6 V6 V6 
(A.8) UE À i > 3 eG ean 
V2, VS 
i 1 
pra DE = eto 0 
2 2 V2 
È A la 0 0 0 


V3 SIN See 


This is the transformation matrix between dii T and Ta 
Similarly, from (A.2) it follows 


(A.9) tr (Ti T;) = Ove . 
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However, 7; and 7,, will not be orthogonal one to another because Ty 
contains T; (through 7 = p:T). By means of another choice of T, they can 
be made orthogonal, but in this way, either (A.6) or (A.7) or (A.9) would 
be violated. Nevertheless, this arbitrariness in choosing 7, is not embarrassing: 
‘ being a linear combination of elements of the product representation DÌ x DI, 
À will not contain elements of the D° representation so that terms con- 
taining T will cancel in all our expressions. 

This is also the reason which enables us to take 


(A.10) (TT) = 0 


’ 


N 
in our calculations. This relation is correct for all angular operators A because 
all non-orthogonal terms between 7; and 7,, are due to 7. 


RIASSUNTO (*) 


Si fa una analisi parziale in onde per la produzione di coppie di bosoni su nucleoni. 
Gli operatori angolari corrispondenti, che caratterizzano completamente la dipendenza 
dallo spin e angolare della matrice S, si esprimono con l’ausilio dei polinomi di Legendre 
e si tabulano. Simultaneamente si dà un metodo generale per calcolare gli operatori 
angolari per processi interessanti più di quattro particelle, casi nei quali il metodo 
diretto conduce a calcoli assai laboriosi. 


(*) Traduzione a cura della Redazione. 
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Magnetic and Crystallographic Analysis by Electron Diffraction. 


S. YAMAGUCHI 


Institute of Physical and Chemical Research - Tokyo 
(ricevuto il 19 Febbraio 1959) 
Summary. — The Lorentz effect observable in electron diffraction is 


utilized for measurement of the magnetization of ferromagnetic sub- 


stances. 


In the present study, the martensitic transformed layer of 


austenitic stainless steel was studied magnetically and crystallograph- 


ically by electron diffraction. 


An electron beam is deflected in a magnetic field. This effect, è.e., the 
Lorentz effect is observable in-the diffraction pattern obtained from a ferro- 


magnetic substance. 


In order to know the conditions existing in the present experiment, a test 
piece, of hard steel with the known remanence was studied by electron dif- 


Fig. 1. — The diffraction patterns of hard 
steel and of gold are superposed. The rings 
are eccentric as a result of the Lorentz 
effect. Wavelength, 0.0292 A. Camera 
length, 495mm. Positive enlarged 2.3 times. 


fraction. The thin edge of a usual 
razor blade (10x10 mm) was em- 
ployed. as test piece. An electron 
beam grazed this sharp edge (mini- 
mum thickness of the truncated 
wedge: about 3000 A) to give rise 
to a diffraction pattern. 

In order to measure the deflec- 
tion of the incident beam caused 
by the magnetic field, a process of 
double exposure with a non-ferro- 
magnetic substance was employed. 
A pattern of gold was beforehand 
photographed, and then that of the 
magnetic specimen was superposed 
on it. The wavelength of the inci- 
dent electrons as well as the posi- 
tion of the photographic plate were 
fixed during this process. Fig. 1 was 
prepared with the hard steel and 
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with a gold foil. In Fig. 1, the diffraction rings of the ferromagnetic 
specimen and those of gold are eccentric as a result of the Lorentz effect. 

Fig. 2 was obtained from two non-ferromagnetic substances, copper and 
gold, in the same process as for 


Fig. 1. All diffraction rings are con- 
centric in Fig. 2, since there is no 
magnetic effect on the electrons. 
Fig. 2 corresponds to the successful 
result of a blank test for the ma- 


gnetic analysis. 


The deflection AZ of the incident 
beam is measurable from the ring 
eccentricity, as illustrated in Fig. 3, 


that is 


INVASI VARIA 


Under the present experimental ar- 
rangement we have a relation bet- 


(1) 


Fig. 2. — Copper and gold. All rings are 


va concentric. No Lorentz effect. Waveleneth, 
ween AZ and the magnetic induc- 0.0306 À = 


tion B of the specimen, 


u I 


NA CT de oy pa, 
mo h 
0 0 


where e means the electron charge (1.6 x10-*° emu), L means the camera 
length (495 mm), m means the electron mass, v means the velocity of the 
electrons, À means the wavelength of the electrons, À means Planck’s constant 


N 


ne 


Fig. 3. — Illustration for 
measurement of the de- 
flection of the incidence. 


(6.6 x 10-27 erg s), and / means the distance travelled 
by the electrons in the magnetic field. 

In Fig. 1 (A=0.0292 A and AZ=2.0 mm) we can 
calculate the mean value of / according to Eq. (1), if 
we assume a mean uniform field for B whose intensity 
is equal to the remanence of the hard steel (10000 G). 
Thus we obtain 


(2) 1=5.7-10¢ A. 


This value was applied for the analysis of the re- 
sults obtained in the following experiments. 

The same process as for the hard steel was carried 
out for a rolled plate of nickel. The edge of the blade 
(10 x10 x 0.5 mm) was sharpened by mechanical poli- 
shing (minimum thickness: about 3 000 A). An electron 


beam grazed this edge of the specimen. Fig. 4 is a double figure consisting of the 
diffraction pattern of nickel and that of gold (A=0.0293 A and AZ=0.66 mm). 
From this figure we can calculate the magnetic induction B of the specimen 
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with the {-value of (2). We obtain 


B= 3300 G. 


This value coincides satisfactorily with the known remanence of nickel. 
When a surface of austenitic 18-8 stainless steel is mechanically polished, 
a martensitic ferromagnetic Beilby layer is formed on the surface as the result 


Fig. 4. — Nickel and gold. Wavelength, Fig. 5. — 18-8 steel mechanically poli- 
0.0293 A. shed and gold. Wavelength, 0.0294 A. 


of the strain-induced transformation (1). This magnetic state was studied 
by the present process. The edge of a test piece (10 X 10 X 0.5 mm) was sharp- 
ened by mechanical polishing. The minimum thickness of the wedge was 
about 3000 A. This edge was ob- 
served by electron diffraction in the 
same way as for Figs. 1 and 4. 
Fig. 5 contains the diffraction pat- 
tern of the test piece and that of 
a gold foil. It is noticeable in Fig. 5 
that the diffraction rings from the 
test piece and those of gold are not 
concentric (A=0.0294 A and AZ = 
= 0.43 mm). It is also recognized in 
Fig. 5 that the series of the rings 
from the test piece is characteristic 
of the body-centred lattice (lattice 
constant=2.86 A). 
A test piece of 18-8 stainless 
Fig. 6. — 18-8 steel electro-polished and steel, which contains a small quan- 
gold. Wavelength, 0.0301 A. tity of molybden (about 2% by 


() J. T. BurweLL and J. WuLFr: Trans. Am. Inst. Mining Met. Engrs., 185, 
486 (1939); S. YamaGUuCcHI: Journ. Chem. Phys., 27, 1310 (1957). 


1200 


(i 28 


MAGNETIC AND CRYSTALLOGRAPHIC ANALYSIS By ELECTRON DIFFRACTION 


bo 
co 
ive} 


weight) as retarder for the transformation, was examined by the present 
process. The ring eccentricity here measured (A=0.0288 A and AZ=0.33 mm) 
was smaller than that in Fig. 5. The diffraction pattern of the test piece 
was characteristic of the body-centred lattice. 

Fig. 6 was obtained from the test piece of stainless steel completely finished 
by electro-polish and from a gold foil. There is found a slight eccentricity 
of the diffraction rings here (A=0.0301 A and AZ=0.17 mm). This magnetic 
effect is not due to the Beilby layer, but due to the ferromagnetic inclusions 
originally existing in the substrate. 


TABLE I 
Material | Polish | Lattice B (gauss) 
| = = | La 
| | | 
Hard Steel | as b.c.c. 10 000 
Rolled Nickel | — | f.c.c. | 3300 
18-8 | Mechanical | b.e. | 2 100 
18-8-Mo | Mechanical | b.c. | 1700 | 
18-8 | Electrolytie | f.c.c. | 800 | 
b.c.c: body-centred cubic 
f.c.c: face-centred cubic 
b.c.: body-centred 


According to Eq. (1) and to the {-value in (2) we can estimate the mag- 
netic inductions B of the test pieces. The results obtained are summarized 
in Table I. The B-values depend on the thickness of the Beilby layer, on 
the orientation of the magnetic domains, and on the nature of the ferro- 
magnetic substance. 


RIASSUNTO (*) 


L'effetto Lorentz osservabile nella diffrazione degli elettroni è stato utilizzato per 
misurare la magnetizzazione di sostanze ferromagnetiche. Nel presente lavoro si sono 
studiati per mezzo della diffrazione degli elettroni i caratteri magnetici e cristallografici 
di uno strato di acciaio inossidabile austenitico che aveva subito una trasformazione 
martensitica. ; 


*) Traduzione a cura della Redazione. 
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(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lascrata 
dalla Direzione del periodico ai singoli autori) 


Electrical Conductivity of Metals at Low Temperatures. 


I. SUPEK 


Institute « Rudjer Boskovié » - University of Zagreb 


(ricevuto il 10 Aprile 1959) 


A previous paper (1) has shown that Bloch’s integral-equation (?) for metal- 
electrons in an eletrical field can, at very low temperatures, be reduced to a dif- 
ferential equation on the energy surface E(k,, k,, k,)— &. The same can be done 
for Peierls’ system (*), where the change of lattice quanta distribution is taken 
into account. 

In an external electrical field there are two stationarity conditions for the; 
electron and the lattice quanta distribution 


of of ON 
(1) | = Ox -—— UE 
ot collisio 8 ot field ot collisioas 
We have applied the same functions for electrons and phonons as Peierls did: 
dfo 
= for ==D(k 
Hk) = fo dE (Fe) 
(2) ; si ote 
Na) dN 1 rq) ho H—é 
= = 5 X= — , lese ti 
DE A RON HN 505 


Hence the stationarity condition takes the form (4) 


eri 15 
(3) cm SÉ oa a ne DL DS 
+ CO (Er — Ex + ho) : [®(k —q) — D(k)- r@il= 
(e-€ + 1)(e-* + 1)(e7— 1) | 
ef 0E df, 
Ca Rena 


(1) I. SuPEK: Zeits. f. Phys:, 149, 324 (1957). 

(*) F. BLocH: Zeits. f. Phys., 52, 555 (1929). 

(3) R. PEIERLS: Ann. d. Phys., 12, 154 (1932). 

(4) Compare all the calculations with the article by SOMMERFELD — BETHE in Mandbuch der 
Phys. Vol. 24/II or WILSON: Theory of metals. 


1202 


ELECTRICAL CONDUCTIVITY OF METALS AT LOW TEMPERATURES 291 


1 
C?Q(E E h Ss - 
ee) rye uen 


[®(k + q) — Ok) — T(q)] = 0. 


It is useful to introduce polar-coordinates with the axis in the direction 


| of grad Æ in the integral over g-space. We can first integrate over # taking into 


account only the factor Q 


f 1 
A sin ds = = QdE= : ok à 
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At very low temperature it is possible to develop the electron function in small 
quantities g,, 9, and q,, the projections of q in the direction of grad E and two 
perpendicular tangents on the surface. 
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In a similar way also 1/grad £ has to be expanded. 
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The averages over ® can be easily calculated (5). 
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The average value in the equation (3) contains the powers in q?, 43, etc. 
We maintain only the part in-g?. The higher powers are negligible at low 
temperatures. By neglecting the change of lattice quanta distribution in the elec- 
trie field and after the integration over E (the factors df,/d# on the left and right 
side act as d-function with the singularity on the H=¢) we obtain the differential 
equation, which can be transformed in ortogonal surface variables. 
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Here A, is the Laplace operator on the surface (A= 4,+A,). These equations 
can be also written in general surface variables (°) 
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~(*) I. SUPEK: Zeits. f.| Phys., 117, 125 (1941). 
(5) GLASER and JAKSIG: Glasnik mat. fiz. 12, 257 (1957). 
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This is just Bloch’s stationarity condition translated into the differential 
form; it does not contain the change of lattice quanta distribution. We have in 
fact to solve the symultaneous system (3) for f(k) and N(q). Here we can first i 
eliminate the phonon’s function by integrating the second equation over the energy | 
surface H(k)._ = Ri 


df 
(8) lsat [D(k + q) — Ok) — T(q)] = 0. 


The phonon’s function must therefore have the same expansion in q,, 9, and q, as 
the electron’s functions 
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By elimination of / we finally obtain the integro-differential equation; in which 
we have to expand the differential forms 
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The transformation in general surface variables is obvious. 
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